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FOREWORD 


This  program  was  sponsored  by  the 


Air  Force  Rocket  Propulsion  Laboratory 
Director  of  Science  and  Technology 
Air  Force  Systems  Command 
Edwards,  California,  93523 


The  technical  effort  reported  herein  was  accomplished  under  Contract 
f 046 1 1  - 71 -C- 0046  and  covered  the  period  from  May  1971  through  June  1973. 

Mr.  Norman  Walker  was  the  Air  Force  project  engineer  at  the  start  of  this 
contract  until  he  left  AFRPL.  Dr.  Randy  L.  Peeters  was  the  Air  Force  materials 
engineer  for  the  remainder  of  this  contract. 


The  program  was  completely  successful  in  meeting  its  objectives.  The 
success  of  the  program  was  due  to  the  combined  efforts  of  key  personnel  at 
two  facilities,  the  Aerojet  Solid  Propulsion  Company  and  the  Texas  A&M  Research 
Foundation.  The  Aerojet  team  performed  an  extensive  experimental  effort, 
developed  computerized  characterization  techniques,  and  performed  theoretical 
work  in  developing  constitutive  equations  which  accurately  model  the  multi - 
axial  behavior  of  composite  propellants,  me  team  at  the  Texas  A&M  Research 
Foundation  supported  the  erfort  by  performing  theoretical  and  experimental 
work  to  develop  mathematical  models  for  characterizing  solid  propellants 
under  realistic  states  of  strain.  The  emphasis  of  their  work  was  to  develop 
simple  constitutive  equations  based  on  thermodynamic  principles  for  visco¬ 
elastic  materials  having  microstructural  damage. 


The  key  technical  personnel  on  this  program  were:  Dr.  Richard  J.  Farris 
who  directed  the  Aerojet  effort  and  coordinated  the  overall  effort  as  the 
Program  Technical  Manager,  and  Dr.  Richard  A.  Schapery  who  directed  the  Texas 
A&M  effort.  Mr.  Frederick  H.  Davidson  performed  all  of  the  experimental  work 
at  Aerojet  and  Mr.  Dennis  F.  Vronay  was  primarily  responsible  for  the  computerized 
characterization  work.  Mr.  L.  E.  Lewis  and  Mr.  R.  T.  Shankle  performed  the 
Texas  A&M  experimental  effort  which  was  under  the  guidance  of  Mr.  Scott  Beckwith. 


This  technical  report  has  been  reviewer!  and  is  approved. 
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Dr.  Randy  L.  Peeters  (MKPB) 
Materials  Engineer,  AFRPL 
Edwards,  California 
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ABSTRACT 


This  program  was  designed  to  develop  three  dimensional  nonlinear 
viscoelastic  equations  that  could  describe  the  stress-strain  response  of 
solid  propellant  materials  for  complex  loading  conditions.  The  approach 
used  was  to  identify  and  mathematically  model  the  underlying  mechanisms 
contributing  to  the  constitutive  nonlinearity  In  these  highly  filled 
Polymeric  materials  and  to  then  include  these  effects  within  the  frame¬ 
work  of  a  continuum  constitutive  theory.  Also  considerable  theoretical 
work  was  done  showing  how  the  resulting  mathematical  representations 
could  be  contained  within  the  framework  of  a  thermodynamic  theory  and 
viscoelastic  fracture  mechanics.  During  the  course  of  this  contract 
Considerable  experimental  as  well  as  theoretical  work  were  performed 
since  the  development  of  the  theory  was  based  on  modeling  observable 
effects.  The  experimental  work  consisted  of  determining  the  three  dimen¬ 
sional  strain  response  due  to  complex  stress  histories  at  temperatures 
from  -65°F  to  +150°F  and  superimposed  hydrostatic  stress-states  from  0 
to  1000  psi.  In  total  over  250  experiments  were  run  using  uniaxial  and 
biaxial  volumetric  dilatometers  and  the  stress-strain-strain  invariant¬ 
time  history  of  each  experiment  is  stored  on  magnetic  tape  and  can  be 
made  available  for  others  to  use. 

The  program  was  successful  in  meeting  its  objectives  and  some  of 
the  program  developments  and  conclusions  are  as  follows.  The  dominant 
mechanisms  leading  to  the  nonlinear  stress-strain  response  is  micro- 
structural  damage  in  the  form  of  the  Mullins'  stress-softening  effect 
and  vacuole  dilatation.  Reversible  nonlinearities  such  as  second  order 
hereditary  strain  effects  were  found  to  be  of  secondary* importance. 

These  effects  have  be >n  successfully  modeled  and  combined  into  a  relatively 
simple  constitutive  tleory  that  works  very  well  for  fitting  the  distor- 
tional  stress-strain  behavior  under  complex  loading  conditions.  This 
theory  is  a  permanent  memory  constitutive  theory  and  contains  irreversible 
effects  of  the  past  history  on  the  current  response  which  are  not  included 
in  the  usual  fading  memory  viscoelastic  theories.  Tin'-  -ffects  are  included 
in  the  permanent  memory  theory  through  time-dependent  structural  damage  as 
well  as  the  ordinary  viscous  energy  dissipation,  me  bulk  effects  were  also 
modeled  and  good  agreement  could  be  obtained  only  when  the  vacuole  gas  phase 
compressibility  was  included  in  the  constitutive  theory  which  resulted  in  a 
mixed  stress-invariant-strain  invariant  representation.  To  handle  the  large  masses 
of  data  generated  on  this  contract  computerized  characterization  techniques 
were  developed  wherein  the  equations  could  be  fit  to  large  masses  of  data 
to  determine  the  applicability  of  the  theory.  For  constant  temperature 
conditions  the  theory  could  be  fit  within  a  standard  deviation  of  +  12%  of 
the  observed  distortional  stress  data  for  approximately  forty  complex  experi¬ 
ments  using  the  entire  response  curve  to  failure.  The  overall  accuracy 
dropped  when  all  the  temperature  data  was  analyzed  together  to  about  +  15%. 

The  hulk  response  predictions  were  poorer  and  typical  deviations  were  +  25% 
nevertheless  all  of  the  proper  trends  were  in  the  predictions  and  much  of 
these  errors  were  no  doubt  experimental. 
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FOREWORD  TO  APPENDIX  A 


Appendix  A,  "Studies  on  the  Nonlinear  Viscoelastic  Behavior  of 
Solid  Propellant,"  represents  the  final  technical  report  from  Texas  A&M 
University  subcontract  effort  to  the  Aerojet  Solid  Propulsion  Company. 

The  work  at  Texas  A&M  was  performed  under  the  direction  of  irofessor 
R.  A.  Schapery.  The  intent  of  this  subcontract  was  to  (1)  support  the 
theoretical  and  experimental  efforts  being  conducted  at  Aerojet,  (2)  pro¬ 
vide  an  independent  laboratory  and  technical  staff  to  test  and  verify  con¬ 
stitutive  assumptions,  (3)  study  the  relation  between  Farris'  .non-fading 
memory  constitutive  theory  and  Schapery' s  thermodynamic  theory,  and 
(4)  compare  the  relation,  if  any,  between  Ferric'  model  fcr  the  Mullins' 
effect  and  Knausa'  theory  for  the  influence  of  microstructural  cracks  on 
the  mechanical  response  of  propellants. 

The  work  performed  at  Texas  A&M  not  only  provided  support  for  the 
work  and  representations  used  at  Aerojet  but  recommended  alternate  Inverse 
representations  using  Lp  norms  in  stress  rather  than  strain.  In  particular, 
the  energy  approach  used  by  Schapery  was  instrumental  in  suggesting  simplified 
three  dimensional  representations  using  only  two  sf-^in  invariants,  the 
dilatation  and  the  octahedral  shear  strain  and  bringing  out  the  coupled 
dependency  of  dilatation  and  mean  pressure  on  octahedral  shear  strain. 
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STUDIES  ON  THE  NONLINEAR  VISCOELASTIC 


BEHAVIOR  OF  SOLID  PROPELLANT 

R.  A.  Schapery 
May  1973 

I.  INTRODUCTION  AND  SUMMARY 

The  state-of-the-art  of  nonlinear  viscoelastic  ciiaracterization  of 
solid  propellant  existing  just  prior  to  the  start  of  this  Investigation 
was  reviewed  In  [1].  Although  some  proolsing  models  had  been  proposed, 
such  as  Farris'  representation  of  microstructural  damage  [2],  little  ex¬ 
perimental  verification  existed  for  propellant  with  vacuole  dilatation. 

The  present  study  was  undertaken  In  cooperation  with  Aerojet  Solid  Propul¬ 
sion  Company  in  order  to  develop  mathematical  models  for  characterizing 
solid  propellant  under  realistic  states  of  strain.  Emphasis  of  the 
study  at  Texas  A&M  University  was  on  the  development  of  simple  consti¬ 
tutive  equations  for  propellant  having  very  large  amounts  of  vacuole 
dilatation  and  on  the  application  of  viscoelastic  fracture  mechanics 
to  constitutive  theory. 

In  Section  II,  constitutive  equations  based  on  thermodynamic  princi¬ 
ples  are  examined  for  elastic  and  viscoelastic  media  with  microstructural 
damage.  These  equations  provide  the  framework  for  establishing  some 
special  cases  which  apply  to  solid  propellant.  The  role  of  fracture 


mechanics  in  constitutive  theory  is  brought  out  using  the  theory  developed 
in  13]. 

The  experimental  investigations  described  in  Section  III  cover  uni¬ 
axial,  strip  biaxial,  and  poker  chip  tests.  Dilatation  was  measured  in 
all  three  modes;  a  dilatometer  was  used  for  the  former  two  tests  and 
changes  in  circumference  were  monitored  to  obtain  dilatation  in  the  poker 
chip  test  [4].  The  strip  biaxial  and  uniaxial  data  were  found  to  support 
the  proposed  constitutive  model.  Poker  chip  test  results  turned  out  to 
be  of  rather  limited  quantitative  use,  especially  with  compressive  load¬ 
ing,  because  of  the  strong  nonlinearities  exhibited  by  the  propellant. 


II.  THEORETICAL  INVESTIGATIONS 


Three-dimensional  constitutive  equations  for  nonlinear  media  with 
microstructural  damage  are  developed  in  this  section.  Results  from 
thermodynamic  theory  are  used  to  establish  the  general  form  of  the 
equations,  while  possible  measures  of  damage  are  examined  in  the  light 
of  viscoelastic  fracture  mechanics  theory.  As  a  means  of  clearly  bring¬ 
ing  out  certain  physical  aspects  of  the  problem  we  first  derive  equa¬ 
tions  for  damaged  elastic  media.  Viscoelastic  equations  are  then  de¬ 
veloped.  In  order  to  provide  a  sufficiently  general  theory  for  incor¬ 
porating  physical  models  of  microstructural  damage  and  to  be  able  to 
examine  the  relative  merits  of  different  specific  representations,  we 
examine  both  the  case  in  which  stresses  are  expressed  as  functions  of 
strains  and  in  which  strains  are  functions  of  stresses. 

For  simplicity  of  notation,  infinitesimal  strain  theory  is  used  in 
the  theoretical  development  of  constitutive  equations.  However,  at  the 
end  of  Section  II  we  show  how  the  results  can  be  lssnediately  extended 
to  finite  strain  theory  by  simply  making  a  change  in  the  notation. 

A.  Constitutive  Equations  for  Elact  Media  wit’'  Microstructural  Damage 

Introductory  comments:  Consider  first  a  nonlinear  elastic  body  with 
an  arbitrary  distribution  of  cracks  with  £l/ed  length.  Derivation  of 
the  overall  stress-strain  equations  can  be  accomplished  by  following  the 
same  steps  used  for  uncracked  bodies.  Specifically,  by  combining  the 
First  and  Second  Laws  of  Thermodynamics,  with  strains  and  temperature 
being  the  independent  state  variables,  there  results  for  isothermal 
processes  [  5  ] : 


(1) 


3F 

e 

0..  “  7 - 
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where  F  is  the  Helmholtz  free  energy/initial  volume  ("strain  energy"), 
e 

It  is  a  function  of  the  strains,  ,  temperature,  T,  and  a  set  of 
parameters  (6^,  say)  which  are  needed  to  define  the  geometry  (size, 
shape,  orientation,  location)  of  all  cracks.  The  inverse  strain-stress 
equation  for  isothermal  processes  likewise  can  be  derived  from  thermo¬ 
dynamics  by  taking  stresses  and  temperature  as  the  independent  variables 
[  3): 
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where  -F  is  the  Gibbs  free  energy  per  unit  initial  volume  ("complemen¬ 
tary  strain  energy").  The  energy  Fg  is,  in  general,  a  function  of  the 
stresses,  ,  temperature,  T,  and  the  crack  parameters  3^. 

Berry  [  6  ]  has  used  the  one-dimensional  form  of  Eq.  (2)  to  derive 
the  overall  Young's  modulus,  Ee>  of  a  large  linear  elastic  sheet  under 
uniformly  applied  tensile  stress,  ,  and  with  one  crack  perpendicular 
to  the  external  load.  He  obtained: 


-  .3 
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(3a) 


where  E  is  the  uncracked  Young's  modulus,  2a  is  the  constant  crack 
length,  and  A  -  sheet  length  x  width.  Thus, 


E  =  - 
e  e 
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If  there  are  N  noninteracting  cracks  of  length  2a^,  instead  of  one  crack, 

N 

then  a*  in  Eq.  (3b)  is  simply  replaced  by  j  ah.  [  7  ) ;  viz. 

K-l 
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Ee  -  a  +  f-  I  *£)-1  E  (3c) 
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Equations  (1)  -  (3)  are  obviously  valia  for  moving  cracks  in  an 
elastic  body  as  long  as  the  velocities  are  small  enough  for  kinetic 
energy  to  be  negligible.  With  th^s  case  one  must  still  perform  the 
differentiation  in  Eqa.  (1)  and  (2)  with  the  parameters  (e.g.  a^ 


held  constant. 


Predictions  of  time-dependent  crack  lengths  and  more  realistic 
overall  properties  will  be  taken  up  it.  a  later  subsection.  For  now  we 
shall  assume  the  instantaneous  values  of  8  aie  known. 

K 

Equations  for  isotropic  media  with  strains  as  the  independent 
variables:  Let  us  return  to  Eq.  (1)  and  assume  that  it  is  valid  regard¬ 
less  of  the  types  of  damage  (e.g.,  submicroscopic  polymer  ch.iin  failure, 
binder  cracks,  and  dewetting)  which  are  characterized  by  the  parameters 
8^.  For  notatlonal  simplicity,  explicit  dependence  of  free  energy  on 

T  and  0„  will  not  be  shown:  but  this  dependence  is  to  be  assumed  unless 
l\  ~ 

stated  otherwise.  Also,  as  noted  previously,  the  partial  derivatives 
with  respect  to  strain  are  taken  while  holding  damage  (viz.  0^)  constant. 
Furthermore,  we  will  assume  that  the  body,  with  or  without  damage,  is 
statistically  homogeneous  and  isotropic.  (The  body  is  defined  to  be 
statistically  isotropic  if,  after  the  external  loads  are  removed,  there 
are  no  preferred  directions;  the  response  to  a  new  set  of  loads  will  be 
independent  of  the  material's  orientation.  If,  for  example,  all  of  the 
cracks  in  a  body  are  parallel,  it  is  not  statistically  isotropic;  if, 
however,  cracks  and  other  flaws  are  sufficiently  random  in  their  orien¬ 
tation  isotropy  can  be  assumed.) 
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Under  the  above  aa sumptions,  the  free  energy  will  depend  on  strain 
at  roost  through  any  three  Independent  Invariants  of  the  strain  tensor 
If,  for  example,  we  use  the  three  Invariants 

1  I  ' 
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1 3  ~  f 1 J  FJ®  E®i 


(4) 


where  repeated  Indices  Imply  summation  over  the  range  1,  2,  3,  then  Eqs. 

(1)  and  (4),  together  with  the  use  of  the  chain  rule  and  thi  assumption 

F  -  F  (1  ,  I  ,  I  ) ,  yield  (where  6  .  .  is  the  Kronecker  delta) , 
e  e  1  2  3 
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The  damage  parameters  0^,  which  are  Implicitly  included  in  F^,  may 

depend  on  the  history  of  the  invariants  (and  their  current  values) ,  but 

not  directly  on  the  Individual  components  of  the  strain  tensor. 

For  physical  reasons,  it  is  often  desirable  to  use  instead  the 

invariant  I  and  the  two  strain-devlator  Invariants i 
1 

*2  ~  *iJ  *iJ 


1  =  e  , 

3  ij  Jm  mi 


where  e^  are  the  devlatorie  strains, 
eij  5  eij  “  I  6ij  \ 

By  substituting  Eq.  (7)  into  (6)  we  find 


I  -  I  -  i  I2 
2  2^1 

I  -  I  -  II  +  f  I3 
3  3  1  2  9  i 


(7) 
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Equation  (5)  can  be  transformed  to  th<*  deviator  strains  by  using 
definition  (7)  and  by  applying  the  chain  rule  to  the  derivatives  in  (5); 
with  F 


F  (I  ,  I  ,  I  )  we  find 
•  1  2  3 


3F  „ 


3F 


’ij 


-  — -  I  6,.  ^  f  --  e,,  +  3 


ij 


31 


IJ 


3F 


31 


e.  e  . 
im  mj 


(9) 


and 
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where  s  are  the  deviatoric  stresses. 
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and  6  is  the  dilatational  stress, 
0  =  °il 


(10) 
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(12) 


Note  that  when  we  set  i  ■  j  in  Eq.  (9)  and  sum  the  three  equations 
(J  -  1,  2,  3),  the  resulting  equation  is  identically  zero.  There¬ 
fore,  Eq.  (9)  defines  five,  rather  than  six,  independent  equations. 

These  five  equations,  together  with  the  dilatational  Eq.  (.10)  ,  form  a 
set  which  is  equivalent  to  the  set  of  six  equations  defined  by  Eq.  (5). 
However,  they  are  now  separated  into  relations  which  express  distortion, 
Eq.  (9),  and  dilatation,  Eq.  (10). 

For  many  metals  and  polymers  (with  or  without  filler  particles)  ,  it 
is  often  possible  to  ne0j.ect  the  dependence  of  mechanical  properties  on 


I  [  8  ] ;  thlf.  simplification  will  be  used  in  all  of  the  following  dis¬ 
cussion  as  present  indications  are  that  it  is  valid  for  solid  propellant. 

Thus,  from  Eq.  (9)  with  F  -  F  (I  ,  I  ): 
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while  Che  form  of  dllatational  Eq.  (10)  Is  unchanged. 

An  important  point  is  that  if  dependence  on  I3  1b  not  neglected, 
the  deviatoric  stresses  in  Eq.  (9)  are  linear  with  respect  to  the  square 
of  the  strain  tensor,  i.e.  e^  e^j  .  This  dependence  implies  the  deviatoric 
stress  tensor  is  definitely  not  a  homogeneous  function  to  degree  one.  On 
the  other  hand,  without  I3,  Eq.  (13)  shows  that  it  is  at  least  possible 
(but  not  necessary)  to  select  F£  such  that  th;  stress  tencor  satisfies 
this  homogeneity  condition. 

An  elastic  free  energy  function  for  propellant:  In  this  subsection 
we  shall  propose  a  specific  free  energy  function  which  will  be  shown  to 
predict  actual  nonlinear  behavior  of  dewetted  propellant.  Specifically, 
we  assume 


Vld,T  +  /  + 

o  o 


where 


IT  5  Ij  -  3aAT 


and  the  term  3aAT  is  the  volumetric  thermal  expansion  due  to  temperature 
change  AT.  Also, 

k,  -  W  ,  g2  -  <yi2> 

N  »  N  (I_)  ,  N  -  N  (I  ) 

1  1  1  2  2  2 

are  four  nonlinear  material  property  functions,  each  of  which  depends  on 
only  one  invariant. 

Substitute  Eq.  (14)  into  (1C)  and  (13);  thus 


(15a) 
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In  the  linear  range  of  behavior,  these  aquations  must  reduce  to 


7  ■  K  I 
3  e  T 


°ij  "  2  Geeij  (16b) 

where  K e  and  are  the  linear  elastic  bulk  and  enear  moduli,  respec¬ 
tively.  Hence,  K.^  and  G^  can  be  interpreted  as  the  nonlinear  generali¬ 
zations  of  and  Ge>  while  Nj  and  N2  produce  coupling  between  dilata- 
tional  and  shear  phenomena.  The  product  must  be  selected  such  that, 

at  small  strains,  it  is  at  least  of  third  order  in  the  strains  in  order 
that  the  proper  linear  limiting  case,  Eq.  (16), is  obtained. 

Let  us  now  go  one  step  further  and  assume 


K  -  K 
1  e 


N,  •  K  I 
1  e  T 


G,  ■  G 
2  e 
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for  12  <  b2 


N2  "  _a(,/l2  "  b)  for  *2  -  b2 

where  a  and  b  are  non-negative  quantities.  Substitution  of  Eqn. 
(17)  and  (18)  into  (15)  yields,  for  I  >_  b2: 


*If  b  “  0,  then  the  free  energy  function  will  be  positive  definite 
if  and  only  if 


a  <  '2G_/K 
e  e 
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while  the  linear  relations  (16)  apply  for  I?  <  b2. 

It  Is  of  Interest  to  examine  Eq.  (19)  in  some  detail  as  these  rela¬ 
tions  turn  out  to  agree  with  a  large  amount  of  propellant  data.  First, 
observe  that  if  Ge,  Ke>  and  a  depend  at  most  on  constant  damage  para¬ 
meters  6  and/or  are  homogeneous  to  degree  zero  in  the  invariants 

iv 

(through  fj^)  ,  then  the  deviatoric  constitutive  Eq.  (19b)  is  homogeneous 

to  degree  one,  and,  in  this  regard,  it  agrees  with  Farris'  observations 

[9].  However,  the  dilatational  Eq.  (19a)  is  not  homogeneous  to  degree 

one  unless  (i)  I  <  b2  and  Kg  depends  on  ^  as  noted  above  or  (ii) 

1  >>  b2  and  K  and  a  both  depend  on  E„  as  noted  above. 

2  e  & 

If  the  dilatational  stress  0  is  either  known  or  is  negligible, 

Eq.  (19a)  enables  dilatation  to  be  easily  predicted;  viz.,  for  1^  >  b2: 

IT  5  -  3aAT  »  -|r-  +  a(/[J  -  b)  (20) 

The  significance  of  this  result  in  predicting  dilatation  in  three 

common  tests  will  be  discussed  in  Section  II-D. 

Consider  the  simple  shear  test  (e  ^  0)  in  which  AT  “  0  and  the 

12 

propellant  is  constrained  such  that  all  normal  strains  vanish,  which 
implies  IT  =  0.  Equation  (19a)  yields,  for  /jJ  =  /2  |e  |>_b, 

f  -  -Ke  a(/2  |eizj-  b)  .  (21) 
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Also,  Eq.  (19b)  yields  a  *•  a  »  a  ,  and  predicts  that  the  effective 

11  22  33 

shear  modulus  is  Ge>  Thus,  the  shear  strain  produces  compressive  normal 
stresses,  which  is  consistent  with  experimental  results  [10]  •  Moreover, 
if  the  propellant  is  permitted  to  expand  during  shearing,  Eq.  (19b)  pre¬ 
dicts  that  the  effective  shear  modulus  will  be  less  than  that  for  no 
dilatation;  this  result,  again,  Is  consistent  with  propellant  data  [10J. 

Let  us  now  return  to  Eq.  (19b)  and  specialize  it  to  uniaxial  loading. 
Assuming  e  >>|  I^|  ,  we  find 


a  ■  3  G  e  [1-6  — ] 
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where 


aK 


6  £ 


/6  G 


As  long  as 


~eV£n 

1-0  —  3  e 
E11 


(22) 


(23) 


(24) 


predicted  stress  (22)  is  the  same  as  that  reported  by  Farris  [11],  where 
6  “  2  -  3. 


We  have  assumed  the  quantities  Kfi,  G  ,  a,  .aid  b  are  constant  Insofar 
as  the  comparison  with  experimental  data  wa-  concerned.  The  influence  of 
damage  parameters  ^  on  these  quantities  will  be  discussed  in  Section  II-E. 

It  should  be  pointed  out  that  the  relations  in  Eq.  (15),  instead  of 
those  in  Eq.  (19),  may  be  needed  to  characterize  propellant  subjected  to 
moderate  and  high  pressures.  Alternatively,  use  of  the  method  described 
in  Section  II-F  together  with  Eq.  (19)  may  prove  to  be  better. 
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Equations  for  isotropic  media  with  stresses  as  the  Independent 
variables:  Constitutive  equations  corresponding  to  the  isotropic 
version  of  Eq.  (2)  are  completely  analogous  to  those  given  in 
terms  of  strains.  We  introduce  the  dilatational  stress  invariant 


9  =  a 


ii 


and  the  two  deviatoric  stress  invariants 
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(25b) 
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The  deviatoric  strains  are  found  to  be 
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and  the  dilatational  strain  I  is 
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(26) 


(27) 


If  Fb  is  independent  of  the  third  stress  invariant  J  then  F.q.  (26) 


becomes  simply 
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(28) 


B.  Constitutive  Equations  for  Viscoelastic  Media  with  Microstruc- 
t.ural  Damage 

Equations  for  Isotropic  media  with  strains  as  the  Independent 
variables:  The  above  considerations  can  be  generalized  to  include  visco¬ 
elastic  behavior  by  drawing  upon  the  irreversible  thermodynsaic  theory  in 
[12].  As  before,  we  will  assume  that  the  dependence  of  material  properties 
on  I  can  be  neglected.  Furthermore,  the.  strain  measures  q^.-.^g,  in  [12] 
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will  be 

assumed  as 

f ollowe 

q  -  if  + 

e 

q  “ 

e 

1 

11 

4 

12 

q  -  4>  + 

e 

q  ■ 

e 

(29) 
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where  4> 

-  $(1  ,  i  ) 

is  a 

material  property  which  is  assumed 

to  vanish 

in  the  absence  of  vacuole  dilatation;  at  small  strains,  this  function 
must  be  at  least  of  second  order  in  the  strains  in  order  for  the  measures 
of  strain  to  reduce  to  the  strains  in  the  linear  range  of  behavior.  The 
significance  of  <p  will  be  seen  by  examining  the  resulting  constitutive 
equations. 

The  use  of  Eq.  (29)  in  the  theory  in  [12]  yields,  after  much  re¬ 
arrangement  , 
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where 

q  -  q(IT,  I  )  E  3*  +  iT  (31) 

is  the  nonlinear  dilatational  strain  measure  which  may  be  a  function  of 

both  I,_  and  I  .  Also,  AG  and  AK  in  these  relations  are  the  LVE  transient 
1  2 

components  of  the  shear  and  bulk  relaxation  moduli,  respectively.  In 
general,  there  may  be  four  different  nonlinear  material  property  functions 
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in  Eq.  (30);  vij. ,  F  ,  a„,  q  (or  $) ,  and  a  ,  where  a  appears  in  the 
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definition  of  reduced  time  <|i, 


iKt)  =  f  dt/a  ;  t|>'  2  i|i(t) 
c 


ae  -  ae(IT,  I2,  T)  (32) 


The  term  in  Eq.  (30b)  involving  AK  gives  rise  to  dependence  of 
deviatoric  stress  on  volume-change  history,  even  when  the  deviatoric 
strain  is  constant. 

Thase  same  equations  could  have  been  derived  more  directly  from  the 
underlying  thermodynamic  theory,  Eq.  (19)  in  [12],  without  making  use  of 
the  constitutive  equations  in  [12]  in  terms  of  strains  and  stresses 
This  alternative  approach  would  consist  of  assuming  that 
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and 

q  2  q  «  q(I  ,  I  ) 
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where  all  seven  q^  are  treated  as  independent  state  variables.  The 
virtual  work  dW,  which  is  used  to  derive  the  constitutive  equations, 
is  written  in  the  form, 


6W  -  s^de^  +  ±  e6IT  +  id 


(34) 


It  is  to  be  noted  that  d^e^  =  e^  2  0,  and  therefore  out  of  the  six 
e^,  only  five  are  actually  independent;  however,  all  six  can  be  treated 
as  independent  quantities  in  the  virtual  work  condition  (34)  as  long  as 
the  term  Ade^  is  added  to  the  equation  (where  X  is  a  Lagrange  multi¬ 
plier)  . 

The  rain  point  we  want  to  make  here  is  that  constitutive  equations 
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(30)  can  be  viewed  as  being  based  on  the  assumption  that  the  thermodynamic 


measures  of  devlatoric  strain,  . . .  ,  are  equal  to  the  deviatoric 

1  6 

strains  themselves,  while  the  tnermodynamic  measure  of  dilacation,  q  ,  is 

7 

a  function  of  the  actual  dilatation,  1^,  and  the  deviatoric  strains 

(through  I  ). 

2 

The  present  constitutive  theory  (30)  can  be  simplified  even  further 
by  assuming  that  when  all  strains  are  constant  in  time  (generalized 
relaxation  test)  the  form  of  the  constitutive  nonlinearity  is  independent 
of  time. 

For  constant  strains  applied  at  t  -  0,  and  using  the  relations 


AK(i|<) 

=  K(i|0  -  Ke 

(35a) 

AG  U>) 

3  G(i/i)  -  iie 

(35b) 

where 


Ke  =  K(»)  and  Gg  5  G(»)  , 

we  find  that  Eq.  (30)  reduces  to 
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(36b) 


The  nonlinear  form  of  Eq.  (36)  at  short  times,  ip  =  0,  will  be  different 
from  that  long  time,  «>,  unless  the  sums  in  square  brackets  vanish; 
viz. , 
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Define  a  new  function  g, 


g  =  q2  +  2  ^  i 


and  then  combine  Eqs.  (37)  and  (38)  to  obtain  a  single  equation  for  g: 

o  .  (40) 

3IT  31  31  3IT 

2  2 

This  equation  is  solved  by  a  standard  method  [13]  to  yield  g  -  g(Fe). 
Namely,  g  is  an  (almost)  >i  bitrary  function  of  the  single  quantity,  Fe 
in  the  linear  range  of  behavior  q  must  reduce  to  Tj,,  which  implies 
g  must  be  such  that 

g  -  |  Fe  when  |Fj  «  1  (41) 

e 

which  is  the  only  restriction  cn  the  function  g(Fe),  From  Eq.  (39), 


q2  -  s<Fe>  -2  r  *2 

e  * 


and  from  Eq.  (37), 


Ke<dg/dFe) 


and  Eq.  (36)  for  constant  strains  becomes 
I  K (it) 
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and  Eq.  (30)  for  transient  strains  becomes 


I  "  [  K<*  “  ♦’>  S  dT 


s1  -  2a^[/  G(ij/  -  i|/')  dr 


+  ^4-  e.  /  K(4i  -  \p')  ^  dt] 


where,  in  Eqs.  (43)  and  (44),  a^,  and  q  are  to  be  expressed  in  terms  of 
the  auxiliary  function  g,  which,  in  turn,  is  an  arbitrary  function  of 
Fe  (except  for  condition  (41)). 

Thus,  the  constitutive  equations  new  contain  three  independent  material 
property  functions:  Fe,  g,  and  a£,  each  of  which  may  depend  on  a  set  of 
damage  parameters,  0^. 

As  one  further  simplification,  assume 

K0|,)  -  A'G(<|i)  (45) 

where  A'  is  a  constant;  this  behavior  was  reported  in  [  4  ]  for  propellant, 
as  determined  from  poker-chip  tests.  Relaxation  relations  (43)  become 
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When  Eq.  (44)  is  specialized  to  a  uniaxial  relaxation  test,  and 


condition  (45)  is  used,  we  find  that  the  dilatation  is  independent  of 
tine;  this  behavior  for  one  propellant  was  reported  by  Farris  [14], 
However,  for  the  PBAN  propellant  studied  on  this  program,  as  well  as 
for  an  82  wtX  CTPB  propellant  studied  at  the  Naval  Ordnance  Station  [15], 
some  change  in  volume  occurred  during  relaxation.  Equation  (44)  is 
sufficiently  general  to  predict  such  behavior,  as  long  as  condition  (45) 
ip,  not  invoked. 

Finally,  attention  is  called  to  the  fact  that  Eq.  (44)  is  capable 
of  predicting  a  history-dependent  pressure  build-up  due  to  shear  strain¬ 
ing,  even  when  1^,  ■  0.  This  behavior  is  a  result  of  the  dependence  of 


Equations  for  isotropic  media  with  stresses  as  the  independent 


variables : 


From  the  thee  •  in  [16]  we  deduce  a  set  of  constitutive 


equations  which  are  similar  to  tnose  in  Eq.  (30);  viz. 
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t  d(s  ,/a  ) 

/  AJ  ( -  ip' )  — -  dx 


(47a) 


t  d(Q/a  ) 

/  AB ( — pj - dx 


(47b) 


where  Fg,  Q,  and  a^  are  material  functions  of  the  invariants  9  and 
J2>  Also  AB  and  AJ  are  the  transient  components  of  the  linear  viscoelastic 
creep  compliances  in  bulk  and  simple  shear,  respe'.tiveiy .  Special  forms 
of  Eq.  (47)  analogous  to  those  in  Eq.  (44)  can  be  dedLxed;  but,  we  shall 
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noc  pursue  these  cases  here. 

C.  On  Constitutive  Equations  which  are  Homogeneous  to  Degree  One 

This  subsection  is  concerned  with  the  following  questions:  (1) 
given  constitutive  equations  in  which  the  stress  tensor  0  is  expressed 
as  a  homogeneous  functional  to  degree  one  in  the  strain  tensor  ,  and 
assuming  unique  inverse  relations  exist  (i.e.,  strain  tensor  is  a  unique 
functional  of  the  stress  tensor)  Are  they  homogeneous  to  degree  one?  ; 
and  (2)  under  what  conditions  do  inverse  relations  exist? 

We  shall  prove  first  that  the  answer  to  question  (1)  is  "yes"  And 
then  state  a  general  criterion  as  the  answer  to  question  (2) . 

Write  the  given  functional  constitutive  equations  as 


°ij  "  8ij(Gkl) 


(48) 


where 


8ij<ackl)  ‘•*u<ekl>  (49) 

with  a  -  constant  scalar.  We  assume  unique  inverse  relations  exist,  i.e., 

Eij  ■  £ij(akl)  <50> 

By  definition  of  f  ^  , 


°ij  "  8ij[fkl(amn)] 


(51) 


It  will  now  be  proved  that  f^  is  homogeneous  to  degree  one  in  the 


stresses;  viz. 


fij(cokl)  “  Ctij  ^°kl) 


(52) 


where  c  -  constant  scalar.  Write  Eq.  (51)  for  stresses  o 
°ij '  “  8ijtfkl<amn,)] 


ij 


(53a) 
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iJij  ■  coij 


Sub  cltutu  f 53b )  Into  (53a), 


°ij  ‘  c  Sij[fkl(c°mn)J 


Upon  using  (49),  with  a  -  1/c,  we  find 


°ij  “  8ij[7fkl(camn)]  *  (55: 

In  view  of  equations  (51)  and  (55),  and  the  uniqueneas  assumption 
(i.e.,  for  a  given  stress  tensor,  there  exists  only  one  strain  tensor), 
we  observe, 


—  f  (ca  )  -  f ,  .  (a  ) 
c  kl  mn  kl  mn 


■  ‘£„i  <”«„>  <56b> 

which,  according  to  equation  (52),  completes  the  proof. 

Now,  consider  the  question  (2).  If  the  quantities  g^ ,  Eq.  (48), 
ere  algebraic  we  can  use  the  well-known  result  from  the  theory  of  coordi¬ 
nate  transformations  that  the  inverse  relations  exist  and  are  unique  if  [5] 

i.  The  six  functions  g  are  single-valued ,  continuous,  and 
possess  continuous  J  first  partial  derivatives  throughout 
the  body  of  interest. 

ii.  The  Jacobian  determinant, 

J  =  1 )  »  d°es  not  vanish  at  any  point  in  the  body. 

A  simple  example  for  which  a  unique  inverse  does  not  exist  is  the 
one-dimensional  uniaxial  case  for  which  a  -  a(c)  has  a  maximum  point  at 
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e  -  e^,  say;  for  all  sCralna  e  >_  e^,  a  unique  single-valued  relation 
e  -  e(o)  does  not  exist. 


In  the  actual  case  of  interest,  represents  a  set  of  six 


fu".tlonals;  i.e.,  it  depends  on  strain  history.  However,  this  set  of 
functionals  can  be  reduced  to  algebraic  functions,  which  can  then  be 
subjected  to  the  above  stated  conditions  for  invertability .  This  re¬ 
duction  is  easily  accomplished  (in  principle)  by  approximating  the  strain 
histories  by  piecewise  linear  or  step  functions  of  time.  In  this  form 


the  original  functionals  g^  become  algebraic  functions  of  the  constants 


defining  the  approximate  strain  histories.  If,  for  example,  we  consider 
all  six  strains  and  let  each  strain  history  be  represented  by  ten  con¬ 
stants,  there  will  be  60  x  60  terms  in  the  Jacobian  determinant. 

In  view  of  the  last  comment,  a  completely  general  Btudy  of  inverta¬ 
bility  would  be  impractical.  However,  significant  insight  into  the 
problem  could  be  gained  by  examining  simple  situations,  such  as  the 
stress-strain  equations  defining  the  behavior  of  a  bar  under  confining 
pressure  and  uniaxial  tension. 


D.  The  Dilatational  Constitutive  Equation 

Characterization  of  the  effect  of  vacuole  formation  and  growth  on 
the  relation  between  dilatational  stress,  6,  and  dilatational  strain, 
Ip,  is  an  important  and  very  difficult  part  of  the  overall  nonlinear 
viscoelastic  characterization  problem.  In  this  section,  the  so-called 
dilatational  constitutive  equation  is  therefore  singled-out  for  special 
study.  Recall  that  some  forms  of  this  equation  for  the  elastic  case, 
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e.g.  Eqa.  (10)  and  (19a),  and  the  viscoelastic  case,  e.g.  Eqs.  (30a)  and 
(47 a),  have  already  been  given. 

We  first  define  the  vacuole  dilatation  I  by  the  following  equation: 

lv  =  ri  3aAT  "  3K~  XT  “  3k 

e  e 

•/here  Kg  -  bulk  modulus  in  the  unstressed  state;  voids  may  exist  in  this 
state.  (Note  that,  in  the  linear  elastic  range  of  ber.avior 

Ij  -  3aAT  +  (58) 

e 

and  therefore  1^  -  0.)  Also  by  definition,  I  is  never  negative. 

In  general,  I  can  be  expressed  as  a  function  of  the  history  of  the 
three  stress  invariants  and  temperature  or  the  three  strain  invariants 
and  temperature.  We  chose  the  latter  four  scalars  here,  and  write 

Iv  -  F  {IT  ,  I2  ,  I3  ,  T}  (59) 

Now,  substitute  Eq.  (57)  for  I  and  find 

4  •  lr +  F  ]Ii  ,  l2  ■  '3  •  T)  •  <6°> 

e 

In  order  to  obtain  a  form  which  is  useful  for  application  to  experi¬ 
mental  data  we  solve  Eq.  (60)  for  mechanical  dilatation: 

XT  “  G  4"  *  ^  *  *3  ’  T}  (61) 

e 

where  G  {}  is  a  function  of  the  history  of  its  arguments. 

Let  us  now  examine  some  special  cases  of  Eq.  (61).  Assume  first 
vacuole  formation  and  growth  depends  only  upon  1^,  T,  and  the  maximum 
shear  strain  eg  (instead  of  1^  and  1^  separately);  the  maximum  shear 
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strain  Is  given  by 


1  c  3 

es  ■  -S-1  <62> 

where  e ^  and  are,  by  definition,  the  algebraically  largest  and  smal¬ 
lest  principal  strains,  respectively;  note  that  tg  >_ 0.  Equation  (61)  Becomes 

h  ■ G  {hr  •  es  •  T)  •  (63) 

e 

A  special  case  of  G,  which  turns  out  to  fit  the  general  form  of 
Farris'  uniaxial,  constant  strain-rate  data  [11]  1b  the  piecewise  linear 
algebraic  function 


:  <  b 

s  s 


vT  +  a«(e«  -  b„>  !  e«  1  b»  i  (as*  ha5  0 


where  a  and  b  may  be  functions  of  6,  T,  and  e  .  It  is  seen  that 

8  S  8 

the  material  is  Idealized  such  that  vacuoles  do  not  form  until  the  maxi¬ 
mum  shear  strain,  eg,  reaches  bg;  thus,  large  positive  values  of  6  are  precluded. 

Alternatively,  we  could,  as  in  Section  A,  suppose  that  the  constitu¬ 
tive  function  is  independent  of  I3>  For  easy  reference  we  record  here 
the  special  result  derived  earlier  (see  Eq.  (20)): 


3iT  5  12<b2 


+  a(/f—  -  b)  :  I  >  b2 


Equations  (63)  and  (64)  will  be  called  the  maximum  shear  strain 
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theory  (MSS)  while  Eq.  (61),  without  dependence  on  13>  and  Eq.  (65)  will 

be  called  the  octahedral  shear  strain  theory  (OSS)  because  this  strain  in¬ 
variant  is  proportional  to  We  now  examine  the  implications  of  Eq.  (64) 

for  the  uniaxial,  strip-biaxial ,  and  doubla-lap  simple  shear  tests. 

Uniaxial  test  (MSS);  For  this  test,  with  or  without  a  confining 
pressure,  the  principal  strains  e  ,  and  e3  are 

e  =  eu  m  applied  axial  strain  (66a) 


e  -  e 
2  3 


The  invariants  I,  and  e  are 
1  s 

ri  “  eu  +  2e3  and  es  "  <eu  *  e3)/2‘ 

Substituting  Eq.  (66)  into  (64)  and  solving  for  e3  yields,  for 

es  ±-bs: 

3aAT  +  -  a8bs  -  (1  ~)  ^ 

e 

c  -  - — 

3  (2  +  ag/2) 

Also  for  e  >  b  ,  the  dilatation  becomes 
8  —  s' 


\  ”  3aAT  +  3K  +  4+a  ^eu  "  -  “  aAT  ~  9K  ^ 


When  e  <  b  ,  set  a  -  0  to  obtain  the  dilatation, 
s  s’  s 

Note  that  the  last  term  in  Eq.  (69)  is  the  vacuole  dilatation;  viz, 


T  r  ^8  AT  0, 

*v  "  4+a  [Eu  "  3  aAT  9K  5 

s  e 
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which  applies  when  Che  condition 


c 


8 


>  b 
—  8 


(71a) 


is  satisfied.  Upon  substitution  of  e^, 


e 


u 


+  aAT  + 


8 

9K 


(68),  into  Eq.  (,71a)  we  find 


(71b) 


which  is  simply  the  requirement  that  the  vacuole  dilatation  be  non-nega¬ 
tive. 

Strip-biaxial  test  (MSS):  The  strip  is  assumed  to  be  clamped  in 
place  at  the  reference  temperature,  for  which  AT  -  0.  Further,  the  grips 
are  assumed  to  be  so  rigid  that  the  strain  in  the  length  direction, 
e2,  is  essentially  zero  for  all  loading  conditions.  The  principal  strains 
are,  therefore 


ex  =  eb  -  applied  strain 

-  0 


e3  «  thickness  strain 


The  Invariants  are 


\  ‘  £b  +  e3 


£b  “  £3 


The  basic  dilatational  Eq.  (64)  yields  for  eg  >_  bg: 

5aAT  +  -  a8bg  -  (1  -  ^)  ^ 

e 

c  - - 

3  (1  +  a8/2) 


and 


1 ,  -  +  !r +  Sr  “b  - b.  - !  *“  -  It1 

1  e  8  e 
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(72a) 

(72b) 

(72c) 

(73) 


(74) 


(75) 


where  Che  last  term  Is  the  vacuole  dilatation: 


Xv  “  2  +  a  ^eb  "  bs  “  2  aAT  “  6K  ^ 
s  e 


(76) 


which  is  valid  when  e  >  b  ;  i.e. 

s  —  s 


%  -  bs +  \ aAT  +  It 


(77) 


When  Eq.  (77)  is  not  satisfied,  the  vacuole  dilatation  is  given  by  Eq.  (76) 
with  afl  -  0. 

Simple  shear  test  (MSS):  The  specimen  is  assumed  to  be  clamped  in 
position  at  the  reference  temperature  for  which  AT  -  0.  It  is  also 
assumed  that  the  outer  two  plates  are  rigidly  fastened  together  and  the 
specimen  is  thin  with  respect  to  the  dimensions  in  the  plane  of  shearing 
and  is  long  in  the  direction  of  shearing.  Under  these  assumptions,  the 
specimen  is  in  a  state  of  simple  shear  in  the  plane  of  shearing.  The 
principal  strains  are 


Y/2 

(78a) 

thickness  strain 

(78b) 

-e  -  -y/2 

1 

(78c) 

where 

y  »  applied  shear  strain 

-  shear  displacement/specimen  dimension  normal  to  displacement. 

The  invariants  I,  and  c  become 
1  s 

I  -  e  and  e  -  y/2  .  (79) 

1  2  3 


A- 26 


rnymmmmmrnmmmmmmmmmmmm 


We  find  the  thickness  strain  (for  eg  >_  bg)  from  Eq.  (64): 

e2 ' 3aLT  +lr +  aa(2  •  b8)  (80> 

e 

and,  in  terms  of  the  maximum  principal  strain,  e^: 

I  -  3a AT  +  It-  +  a  (e  -  b  )  .  (31) 

1  3K  sis 


where  the  vacuole  dilatation  Is  now 

I  «  a  (e  -  b  )  (82) 

v  s  i  s 

The  condition  e  >  b  is  the  same  as  e  >  b  :  also,  set  a  t  0  whenever 

c.  —  a  .  —  a  ’  ’  a 


Representation  of  vacuole  dilatation  In  terms  of  strains  due  to 
externally  applied  loads  (MSS) :  If  after  clamping,  but  prior  to  loading 
in  a  test  machine,  the  specimens  are  subjected  to  a  temperature  change,  A1 , 
and  pressurization,  p  »  -0/3,  the  strains  in  this  state  ire: 


e  =*  aAT  +  4m—  ■*  aAT  -  4m—  (uniaxial  test)  (83a) 

U  “K  JtN. 

e  e 

c^  **  aAT  +  aAT  -  (biaxial  test)  (83b) 

e  e 

e  »  0  (shear  test)  (83c) 

1 

where,  for  the  biaxial  and  shear  tests,  we  have  assumed  v  -  which 
implies  the  three  principal  stresses  under  external  pressure,  p,  are  all 
approximately  equal  to  -p. 

If  the  specimens  are  now  strained  in  the  test  machine,  the  principal 
strains  eue>  e^e,  ege>  say,  due  to  the  externally  applied  load  are  equal 
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to  the  difference  between  e^,  and  e^,  and  and  the  respective  initial 
•trains  (83) : 


c  -  e  -  a  AT  +  -r— 
ue  u  3K 

e 

(uniaxial  test) 

(84.) 

ebe  ‘  "  2  aAT  +  IT 

e 

(biaxial  test) 

(84b) 

e  ™  e  i 

se  1 

(shear  test) 

(84c) 

Now,  0  for  the  uniaxial  test  is 

0  «  -3p  +  a 
r  u 

(85a) 

and  for  the  biaxial  test 

3a, 

0  =  _3p  +  ___ 

(85b) 

where  au  and  ob  are  the  stresses  applied  by  the  testing  machine  to  the 

uniaxial  and  strip-biaxial  specimens, 

respectively. 

Upon  substituting  Eqs.  (84)  and 

(85)  into  EqB .  (70),  (76), 

and  (82) 

we  find 

3a  .  o 

I  -  tt2-  [  e  -4b-  q£-] 
v  4+a  ue  3  s  9K 

R  e 

(uniaxial  test) 

(86a) 

2a  „  °u 

Xv  “  '  ^Se  “  bs  “  4K^ 

s  e 

(biaxial  test) 

(86b) 

I  ■  a  ,  c.  -  b  ] 
v  s  se  s 

(shear  test) 

(86c) 

Recall  that  K  is  the  bulk  modulus  in  the  unstressed  state, 
e 

Thus, 

the  terns  °u/9Ke  and  a^/4Ke  are  expected  to  be  relatively  small,  and  there¬ 
fore  will  be  neglected  in  the  remainder  of  Section  D.  Analogous  terms  will 
be  neglected  in  the  OSS  theory. 
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We  now  turn  to  the  application  of  Eq.  (65)  to  prediction  of  dilata¬ 
tion  in  the  three  teste.  Enaentially  the  sane  steps  as  in  the  MSS  theory 
are  followed  in  predicting  vacuole  dilatation  from  the  'Job  theory;  con¬ 
sequently,  we  shall  omit  the  details  of  the  derivations,  itie  strains 

e  ,  e,  .  and  e  are  defined  in  the  same  way  in  both  theories, 
ue’  be’  se 

Uniaxial  test  (OSS):  The  vacuole  dilatation  ia  found  to  oc 
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Figure  2.  Comparison  of  Octahedral  and  Maximum  Shear  Strain 
Theories  for  a,,  =  0.5. 


Figure  3.  Comparison  of  Octahedral  and  Maximum  Shear 
Strain  Theories  for  au  =  1. 
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Is  Che  slope  of  Che  dllaCadon  vs,  strain  line  in  uniaxial  tension;  also, 
is  Che  value  of  uniaxial  sCrain  ac  which  Che  dilatation  Line  inter- 
seccs  Che  sCrain  axis.  The  same  values  of  a  and  c„  wen?  used  in  boch 

Ll  D 

Cheories  so  Chat  Che  same  vacuole  dilacacion  would  be  predicted  for  the 
uniaxial  test.  With  the  normalized  variables  used  in  Figs.  1-3,  the 
graphs  are  independent  of  and  the  predictions  in  Fig.  I  are  indepen¬ 
dent  of  au  as  well;  thi  range  0  <  <_  1  in  these  figures  brackets  the 

range  of  values  a^  observed  for  solid  propellant  under  constant  strain 
rate. 

It  is  seen  that,  for  given  and  e^,  the  following  inequality  is 
true  for  each  theory: 

Iy  (uniaxial)  <  (biaxial)  <_  (shear)  (90) 

Moreover,  except  for  the  shear  prediction  in  Fig.  3,  the  MSS  theory  pre¬ 
dicts  a  greater  dilatation  than  the  OSS  theory  in  each  test.  These  ob¬ 
servations  provide  a  basis  for  evaluating  the  relative  accuracy  of  each 
theory  from  experimental  data.  Note  that,  according  to  these  figures, 
use  of  shear  and  uniaxial  data  is  better  for  this  purpose  than  biaxial 
and  uniaxial  data  when  au  is  close  to  unity. 

E.  Microstructural  Damage  Theory 

In  this  section  the  damage  parameters  and  their  effect  on  constitutive 
equations  are  examined  in  the  context  of  viscoelastic  fracture  mechanics 
theory.  First,  however,  we  shall  consider  their  role  in  the  special  dila- 
tational  constitutive  equations  (64)  and  (65),  and  in  the  shear  equation 
(19b)  with  respect  to  vacuole  dilatation. 
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Vacuole  dilatation:  According  to  Farris'  micros tructural  model  for 
vacuole  dilatation  [11],  the  volume  fraction  of  solids  about  which  vacuoles 
exist  is  proportional  to  the  slope  of  the  vacuole  dilatation  vs.  uniaxial 
strain  curve;  experimental  studies  showed  that  this  constant  of  proportion¬ 
ality  is  approximately  unity  [11].  Thus,  the  Increase  of  dilatation  with 
strain  when  dlv/deu  ■  constant  is  attributed  to  expansion  of  the  cavities 
around  the  particles.  With  this  model  as  motivation,  let  us  now  suppose 
for  the  sake  of  argument  that  (i)  the  relative  change  in  damage  parameters 
is  small  when  dly/de^  is  constant  and  (ii)  this  change  per  cycle  remains 
small  in  subsequent  unloading  and  loading  cycles.  We  are  assuming,  in 
effect,  that  the  major  change  in  B„  occurs  during  the  first  loading  as  the 

invariant  I  (or  e  )  is  increased  from  zero  to  its  value  when  dl  /dc 
2  8  v  U 

first  becomes  constant,  and  that  there  is  little,  if  any,  rehealing  when 
the  speciment  is  unloaded. 

In  Farris'  model  [11]  and  in  the  model  proposed  above  it  is  not 
actually  necessary  to  assume  the  microstructural  damage  is  entirely  in 
the  form  of  dewetting.  The  case  in  which  many  (or  all)  particles  have  a 
skin  of  binder  adhering  to  them  is  not  precluded  and,  as  a  result,  the 
term  vacuole  dilatation,  rather  then  dewetting,  is  used  in  this  report. 

Referring  now  to  the  simple  constitutive  relations  Eq.  (19)  we  see 
that  if  the  above  two  suppositions  are  valid  then  the  parameters  Kfi  and 
Gg  are  the  bulk  modulus  and  shear  modulus,  respectively,  in  the  range 
1^  <  b2,  after  the  sample  has  undergone  strains  for  which  12  >  b2.  Some 
change  in  these  parameters  from  cycle-to-cycle  can  be  expected  as  addition¬ 
al  damage  occurs  with  each  cycle.  The  bulk  modulus  Kg  will  be  less  than 
that  for  the  undamaged  propellant  since  the  specimen,  after  being  unloaded, 
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will  have  a  larger  void  volume  fraction  than  existing  in  the  undamaged 
propellant;  this  point  is  supported  by  the  experimental  data  in  Sec¬ 
tion  III. 

The  above  model  can  be  easily  incorporated  into  viscoelastic  con¬ 
stitutive  equations.  For  example,  referring  to  Eq.  (44),  in  which  all 
three  material  property  functions  depend  on  the  elastic  free  energy  F  , 
we  would  simply  use  for  Fg  the  special  form  in  Eq.  (14)  together  w.th 
Eqs.  (17)  and  (18). 

Microcrack  initiation  and  growth:  The  state  of  the  propellant  at 

strains  for  which  dl  / de  is  constant  is  the  end  result  of  an  extensive 
v  u 

amount  of  micro-flaw  growth  which  occurs  within  the  binder  and/or  be¬ 
tween  binder  and  filler  particles.  In  this  subsection  we  model  the 
initial  flaws  as  cracks  and  examine  the  implications  of  viscoelastic 
fracture  mechanics  theory  in  predicting  the  dependence  of  damage  para¬ 
meters  on  applied  loading  history.  Thus,  while  the  preceding  subsection 
dealt  with  the  final  stage  of  damage,  this  subsection  is  concerned  with 
the  Initial  stage.  Through  a  combination  of  the  two  models  a  constitutive 
theory  applicable  to  all  strains  is  proposed. 

In  the  following  theoretical  development  we  assume,  for  simplicity, 
that  the  viscoelastic  body  is  subjected  to  external  uniaxial  tensile 
stress,  o  -  J(t)  ,  and  that  the  cracks  propagate  in  only  the  so-called  open¬ 
ing  mode.  Proposed  generalizations  to  an  arbitrary  state  of  -.pplied  stress 
and  to  other  crack  propagation  modes  are  given  at  the  end  of  this  sub¬ 
section. 


We  suppose  that  there  exists  an  initial  distribution  of  cracks  within 
the  binder  and  possibly  between  binder  and  filler;  the  size  of  each  crack 


is  assumed  co  be  defined  by  a  single  parameter,  2aQ,  representing  the 
length  or  diameter  of  the  crack.  As  the  loading  on  the  body  Is  in¬ 
creased  from  zero  a  given  crack  will  start  to  grow  when  the  stress  in¬ 
tensity  factor*  at  the  tip,  N  ,  exceeds  a  certain  critical  value,  N 

J  o  or 

This  value  is  identical  to  that  for  an  elastic  continuum  whose  Young's 
modulus,  E,  and  Poisson's  ratio,  v,  are  equal  to  the  long-time  (rubbery) 
values  Er  and  vr  respective)/,  for  the  actual  viscoelastic  material  [3  ]; 
thus 


N 


or 


(91a) 


where  r  is  the  fracture  energy.  If  the  stress  intensity  factor  equals  or 
exceeds  the  value  based  on  initial  (glassy)  modulus  and  Poisson' b  ratio 


(91b) 


nen  the  crack  velocity  is  predicted  to  be  very  high,  and  is  limited  only 

iv  wave  action.  We  shall  say  the  material  surrounding  a  given  crack  is 

tailed"  if  N  >  N  ;  the  time  at  which  N  is  first  reached  will  be 
o  -  og’  og 

uxxea  the  failure  time. 

’’he  fracture  theory  in  [  3  j  can  be  used  to  predict  time-dependent 


"The  stress  intensity  factor  is  defined  in  the  singular  stress  dis- 
oution  near  the  crack  tip. 


o  ”  N  /  /x 

y  ° 


/nett 


is  the  stress  normal  to  the  crack  plane  and  x  is  the  distance  ahead 


crack  size  when  N  <  N  <  N  and  failure  time  for  each  crack  in  a 
or  o  og 

propellant  specimen  under  the  following  assumptions:  (i)  the  rubber 
binder  is  linearly  viscoelastic;  (ii)  the  Poisson's  ratio  of  the  binder 
is  constant  (we  shall  assume  v  -  1/2);  (iii)  all  filler  particles  are 
rigid  relative  to  the  binder;  and  (iv)  the  propellant  specimen  is  linear¬ 
ly  viscoelastic  with  crack  sizes  fixed,  and  is  under  a  spacewise  uniform 
temperature. 

Now  the  crack  velocity,  4,  when  <  Nq  <  N  is  governed  by  the 

equation  [ 3  ] 

D(t  ) - —  (92a) 

3w  N „2 


t  -  a/3a 
a 


a  -  2/a  2I  2  (92c) 

o  m  1 

where  existing  information  on  rubber  indicates  that  the  fracture  energy  r 

and  the  stress  distribution  in  the  failure  zone  behind  the  tip,  a  I  , 

m  i 

are  constant.  The  quantity  D(t  )  is  the  creep  compliance  in  uniaxial 

a 

tension  expressed  in  terms  of  the  "effective  time,"  t^;  also,  a  is  the 
length  of  the  failure  zone  behind  the  tip. 

It  is  further  assumed  that  during  most  cf  the  time  required 
for  local  failure,  the  creep  compliance  is  given  by  the  power  law,* 


D(t)  -  Mt/a^)1 


where  D  and  n  are  constants.  Also,  from  the  time  the  body  is  first  loaded, 

1 

t  -  0,  to  the  failure  time  for  a  given  crack  we  assume  the  crack  is  isolated. 


*However,  at  very  short  and  long  times,  D(t)  -*■  E  1  and  D(t)  -*■  E^, 
respectively. 
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That  is,  the  crack  is  assumed  to  be  so  small  that  the  only  geometric 
parameter  affecting  crack  growth  is  its  own  instantaneous  size,  2a;  in 
support  of  this  point  for  cracks  which  are  initially  isolated,  it  is 
shown  in  [  3  ]  that  the  crack  size  at  90  percent  of  the  failure  time  is 
only  10n  x  (2aQ) ;  e.g.,  if  n  *  0.3,  most  of  the  time  required  for  failure 
Is  consumed  while  the  crack  doubles  in  size.  The  stress  intensity  factor 
for  the  i—  Isolated  crack  satisfying  the  assumptions  (i)  -  (iv)  listed 
Just  above  Eq .  (92)  can  be  written  in  the  form  [  3  ] 

N  -  f .  o  (94) 

ol  1  1 

where  2a^  is  the  size  of  the  i—  crack  and  the  coefficient  ^  1b  a  measure 
of  state  of  stress  existing  in  the  binder  in  the  neighborhood  of  the  1— 
crack  relative  to  the  applied  stress;  f^,  which  will  be  called  a  stress 
concentration  factor,  is  independent  of  material  properties  and  depends 
only  on  the  local  particle  geometry  and  spacing.  According  to  the  theory 
in  [3],  Eq.  (94)  is  identical  to  the  stress  intensity  factor  in  an  elastic 
binder;  thus,  in  principle,  elastic  analysis  can  be  used  to  derive  f^  if 
desired. 

We  combine  Eqs.  (92-4)  and  find  the  instantaneous  size  of  the  i— 

■“V, 

crack  at  time  t: 

3(l-n)  ,(2*1) 

t  _ l _ 

- rr~ - T7 —  “  f  .q  /  .  nr  2n  T  2n 

1/n  1/n  i  4n  T  o  I, 

a  .  a.  o  ml 

ol  i 

where  the  curly  bracket  in  the  integrand  will  be  a  function  of  time  if 
chemical  aging  occurs  and/or  if  there  is  rehealing  (which  are  accounted  for 
through  time-dependent  properties).  Also, 


1/n 

o  q 

—  dt  (95a) 
“T 
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q  *2(1+$ 


(95b) 


The  failure  time  of  this  crack  is  denoted  by  t^,  and  is  obtained  by  letting 


a^-*-  <•;  it  can  be  written  in  the  form  of  a  linear  cumulative  damage  rela¬ 


tion: 


1  -  ;  1  r  -' 
o  ci 


(96a) 


,  th 


where  t  .  would  be  the  failure  time  for  the  i —  crack  if  the  specimen  had 
cl 


been  subjected  to  a  timewise  constant  stress,  a,  (applied  at  t  ■  0)  equal 
to  the  actual  stress  at  the  time  t,and  had  constant  material  and  fracture 
properties  equal  to  those  actually  existing  at  the  time  t: 


_  1  _  1  _  1 
tci  -  (A  n)  (f  i"*1)  <aoi  n)(M  n)(a"q)aT 


(96b) 


where  A  is  a  constant  and  is  the  group  of  terms  in  curly  brackets  in  Eq.  (95a) 
at  a  preselected  reference  state,  and  M  is  the  ratio  of  this  group  at  time 
t  to  the  reference  value.  Without  aging  and  rehealing  M  is  unity;  other¬ 
wise  M  *  M(t),  and,  although  not  essential,  we  assume  M  is  the  same  func¬ 
tion  for  all  cracks. 

We  substitute  Eq.  (96b)  into  (96a)  and  solv'-  for  the  stress  concentra¬ 


tion  factor  required  to  produce  failure  at  time  t^: 


Ba 


oi 


1_ 

nq 


(97a) 


•Mq 


where 


B  5  A 


nq 


and  llallMq  is  a  "weighted"  Lebesque  norm, 


’I  1 


5  >. 

1  . 


(97b) 


in  which  C  is  reduced  tine, 

5  =  /  ^  (97d) 

o  *T 

Of  course,  without  aging  and  reheaiing  we  set  M  -  1  and  obtain  the  LebeBque 
norm  itself  in  terms  of  reduced  time.  Note  also  that,  from  Eq.  (95b), 


i _ 1 _ I.  I 

nq  2(n+l)  2  q 


(97e) 


Thus,  for  all  cracks  in  which  f  >_  f^  failure  of  the  surrounding 

material  has  occurred.  Cracks  with  smaller  concentration  factors  have  not 

yet  failed.  When  the  time  t^  is  exceeded  for  s  given  crack,  the  growth 

will  be  very  rapid  until  the  tips  are  arrested  at  filler  particles  and/or 

the  tipo  move  into  a  region  of  low  stress  concentration.  We  shall  assume 

that  the  increase  in  the  1—  crack  size  from  the  time  t^  to  the  time  the 

crack  is  arrested,  tfli>  say  (or,  at  least,  until  its  velocity  reduces 

to  a  relatively  small  value)  is  much  larger  than  2sq,  and  that  the  time 

difference  tfll  -  t^  is  negligible  compared  to  t^.  Let  us  denote  the  crack 

size  at  time  t  .  as  2a  , :  since  a  .  >>  a  .  and  in  view  of  the  comments 

immediately  following  Eq.  (93)  we  conclude  that  the  influence  of  the  i— 

crack  on  overall  mechanical  response  is  felt  approximately  at  the  time 

t^,  and  the  magnitude  of  this  Influence  is  directly  related  to  the  total 

growth.  2a  ,  -  2a  .  a  2a  , .  Interaction  between  cracks  when  t  >  t.  is 
°  ’  ai  oi  ai  - - i - 

allowed. 
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Let  us  now  predict  the  increment  in  specimen  strain  thai  results 
when  the  i—  crack  grows  to  length  2a^  at  the  time  L^.  In  view  of 
the  linearity  assumption  this  increment  can  be  calculated  by  superposition. 
Specifically,  consider  the  normal  stress  distribution  that  acts  across 
an  imaginary  surface,  where  this  surface  is  identical  to  the  fracture 
surface  which  is  later  formed  when  the  crack  grows  to  length  2aa^.  The 
increment  in  specimen  strain  due  to  actual  crack  growth  is  identical  to 
that  resulting  from  a  sudden  superposition  of  a  pressure  distribution  at 
time  t^,  whose  magnitude  is  equal  to  the  normal  stress  distribution 
existing  across  the  imaginary  surface  following  failure.  This  pressure 
is  proportional  to  applied  stress  a.  Linearity  plus  dimensional  analysis 
yield  the  following  strain  increment  for  a  thermorheologically  simple 
material  If  the  applied  stress  is  constant  for  t  >  t^: 

*  S±  a  D<£  -  ?±)  (98a) 

where  £  is  defined  in  Eq.  (97d)  and 

5iH/tl^  (98b) 

1  o  *T 

Also,  D(£)  is  the  uniaxial  creep  compliance  of  thP  propellant  in  terms  of 
reduced  time;  it  should  be  added  that  with  rigid  particles  and  linearity, 
dimensional  analysis  can  be  used  to  show  that  this  creep  compliance  is 
proportional  to  the  insitu  compliance  of  the  binder.  The  coefficient 
is  a  constant. 

Extension  of  Eq.  (98)  to  a  time  varying  stress  following  time  tfli 
is  readily  accomplished  by  superposition  as  long  as  further  change  in 
crack  length  Is  small  compared  to  2aa^.  Specifically,  introduce  the 
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Heaviside  unit-scep  function, 


H(C  -  Ci)  = 


1,  «  >  5, 
0,  5  <  5, 


where  5  =  C(ti>.  Then  the  strain  is 


d[H(S*  -  CJa] 


ei  “  Si  ;  D<£  '  *’>  dT 
0 


A  '  J 

- -  d€’ 


Since  the  Dirac  delta  function,  6,  is  given  by 

d  H<5  -  5  J 

5(5_5i)._ - —  (101) 

Eq.  (100)  can  be  written  explicitly  as 

Ei  *  SiciD<5  '  V  +  st  '  -  «’>  ffr  iV  (102) 

£i 

where  a =  c(?i). 

We  now  turn  to  a  statistical  description  of  the  microcracking  problem, 
and  thereby  derive  a  constitutive  equation  by  summing  over  the  strains 
given  by  Eq.  (102).  First,  define  g^  as 

8l . 

The  stress  concentration  factor  f^  and  initial  flaw  size  are  random  variables 
which,  in  turn,  imply  is  a  random  variable.  Likewise,  is  a  random 
variable.  Dropping  the  subscript  (i)  we  write  the  expression  for  the 


number  of  cracks,  Nc,  having  g-values  between  g  and  g  +  dg,  and  S-values 
between  S  and  S  +  dS  in  the  form 
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n(S,  g)  d  S  d  g 


(104) 


Nc  ‘ 

The  distribution  function  n  -  n(S,  g)  can  be  calculated  from  a  distribution 

function  in  terms  of  the  three  random  variables  a  ,  f.  and  S,  which  will  be 

o’ 

demonstrated  later.  Equation  (104)  will  be  used  in  the  constitutive  theory 
development. 

It  is  to  be  noted  that  g  and  failure  time  are  directly  related  through 
a  deterministic  relation,  which  is  obtained  by  combining  Eqs.  (97a)  and  (103), 


g(Cf) 


«f  ± 

/  Mn  aqdf; 


1 

q 


(105) 


where  is  now  used  as  the  generic  variable  representing  reduced  failure 
time. 


The  strain  due  to  all  microcracking,  e^,  say,  at  reduced  time  £, 
is  obtained  by  multiplying  Eq.  (100)  by  (104),  replacing  by  F,^,  and 
integrating  over  all  possible  S-factors  (-  »  <  S  <  «)  and  all  possible 
g-values  corresponding  to  0  <_  ±  ?.  There  results,  finally, 


where 


and  G 


e 

m 


€ 

f  D(5 


a  /  G(g') 

g’ 


lg 


1 

V  - 

/  Mn  aq  d5 
o 

G(g)  is  the  following  distribution  function: 


d?’ 


(106) 


(107) 


oo 

G  =  /  S  n(S,  g)  dS  (108) 

—  OO 

In  view  of  the  linearity  assumption,  the  total  axial  strain,  e,  is 
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Chat  due  t"  the  externally  applied  stress  acting  on  undamaged  propellant, 

eu  =  /  D«-{')}|r«'  (109) 

0 

plus  the  strain  cm;  viz,  from  Eqs.  (106)  and  (109), 

4  (  00  i 

e  -  /  D(C  -  V)  oil  +  /  G(g)  dg]  [  d£»  (110) 

o  dH  8*  1 

where  g'  Is  given  by  Eq.  (107)  in  terms  of  a  weighted  Lebesque  norm  at  time 

00 

In  view  of  Eq.  (107),  we  see  /  is  a  non-decreasing  function  of  time 

g' 

(without  rehealing). 

A  distribution  function  G(g)  of  particular  interest  is  the  power  law 

k  g'r  ,  g  <  8, 

G  -  (111) 

(°  *  *  >  *n> 

where  g^,  Gj,  and  r  are  constants;  this  function  will  be  used  later. 

Before  examining  some  special  cases  of  Eq.  (110)  we  shall  derive  the 

relation  between  n  -  n(S,  g)  In  Eq.  (104)  and  the  dist-ibution  function 

for  stress  concentration  factor  and  initial  crack  size.  The  distribution 

function  m  *  m(S,  aQ,  f)  is  defined  by  the  condition  that  is  the  number 

of  cracks  having  S-values  between  S  and  S  +  dS,  a  -values  between  a  and 

a  +  da  ,  and  f-values  between  f  and  f  +  df : 
o  o 

N,  -  m(S,  a  ,  f)  dS  da  df  (112) 

ta  o  o 

The  number  of  cracks  Nc  in  Eq.  (104)  is  equal  to  the  integral  of  over 

the  area  in  the  aQ  -  f  plane  bounded  by  the  curves  g  -  constant  and 
g  +  dg  ■  constant,  where  f  and  aQ  satisfy  Eq.  (103).  After  equating  Nc 


Co  the  area  there  results 

_  1_  _  1_ 

o(S,  g)  -  B  /  m(S,  aQ,  BgaQ  nq)  aQ  r*q  daQ  (113) 

o 

Aa  an  example,  if  m  is  a  power  law  in  ft 

m  -  m  1  (S ,  aQ)  f°  (114a) 

where  a  is  a  constant,  we  find  that  n(S,  g)  is  given  by  a  power  law  in  g, 

.  -  —  (1+a) 

n(S,  g)  “  B  “  g“  /  mj  (S,  aQ)  aQ  nq  daQ  (114b) 

o 

By  substituting  this  result  into  Eq.  (10B)  the  function  G(g)  is  found  to 
be  the  power  law  in  Eq.  (Ill)  with  r  »  -a  and  g^  -  »>. 

We  now  examine  two  rpecial  cases  of  Eq.  (110)  in  order  to  help  bring 
out  the  physical  significance  of  this  constitutive  equation. 

Elastic  specimen:  It  is  assumed  that  the  creep  compliance  appearing 
in  Eq.  (110)  is  a  constant,  D,  say;  but  the  compliance  in  the  fracture 
theory  is  still  the  power  law,  Eq.  (93).  These  assumptions  are  not  incon¬ 
sistent.  This  point  follows  from  the  fact  that  the  effective  time  t  in 

a 

fracture  Eq.  (92)  is  typically  several  decades  smaller  than  that  governing 
the  overall  mechanical  response  of  the  binder  [3];  thus,  for  example, 
crack  growth  may  be  controlled  by  viscoelas'.ic  processes  at  the  same  time 
the  binder  response  is  essentially  in  its  long-time  elastic  response  range, 
Equation  (110)  reduces  to 

!CO  k 

1  +  /  G(g)  dg  ? 

8  ’ 


inhere 


-1 


(115) 


Equation  (115)  shows  that  the  microscopic  damage  is  characterized  by  a 
single  function  of  the  weighted  Lebesqua  norm  ll0l!Mq»  88  in  (HO). 
This  result  is  analogous  to  that  obtained  by  Farris  [  2 ] ,  but  the  role 
of  stress  and  strain  is  reversed. 

The  power  law  distribution  Eq.  (Ill)  reduces  this  theory  to 


c  ■ D0  J1  kIMI*,)'"1-  8-1"r’! 


when  til,  and 


e  ■  Do 


1  +  c^ilMl  g,J 


'Mq  - 


(116a) 

(116b) 
-1 


when  r  ■  1.  Solution  (116)  is  valid  for  only  the  range  | 

when  I  ML  <  g  ^  failure  has  not  occured  and  one  must  set  G.  »  0. 

1 1  1 'Mq  m 

If  -  <»,  the  integral  in  Eq.  (115)  converges  only  if  r  >  1;  for  this 
case  Eq.  (116a)  shows  that  the  strain  will  be  approximately  homogeneous 
to  degree  one  in  stress  if  r  is  close  to  > nity. 

Creep  of  a  viscoelastic  specimen:  We  set  0  "  °c  H(£),  where 
oc  -  constant,  and  use  ths  viscoelastic  theory  to  predict  creep  compliance 
with  microcracking,  Dc  =  e/oc<  For  simplicity,  isothermal  behavior  (aT  -  1) 
without  rehealing  and  aging  (M  ■  1)  is  assumed  and  both  power  laws  Eqs.  (93) 
and  (111)  will  be  used  with  gm  -  ®. 

From  Eq.  (107), 


g'  -  g(t')  -  ac  A(t') 


where  we  have  set  “  t'.  Equation  (110)  yields 


_  e 


G.B 


r-1 


D  =  ~  “  D  tn  [1  +  -i-  (o  )r_1  tq  ] 
cal  q  c 

c 


(117) 


(118) 
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where  the  constant  B  la  a  Beta  function, 


B  5  /  (1  -  u)n  (u)e 


T  (r-1)  -1 


For  propellant  n  =>  0.2,  which  Implies  q  “  12,  and  therefore  if  r  la  close 
to  unity  we  see  -  tn  and  that  D£  is  nearly  independent  of  stress.  The 

value  r  ■  1  is  not  admissible  since  B  ■  »  for  this  case. 

Generalizations  of  the  constitutive  theory  with  microcracking:  The 
above  constitutive  theory,  Eq.  (110),  is  a  special  case  of  the  nonlinear 
theory,  Eq.  (47),  based  on  thermodynamics.  Indeed,  we  may  recover  the 


present  theory  from  Eq.  (47)  by  assuming: 


(120a) 


30  “  9^9aG 


(120b) 


3Fd 

—  -  J(-)/4a_ 
G 


(120c) 


(120d) 


—  -  1  +  /  G(g)  dg 
®G  s' 


(120e) 


where 


E*  ±  j-  i 

/  m”  j  AjB  +  A2/J2|q  dd  q 


(120f) 


and  A  (or  A  )  is  an  arbitrary  constant.  However,  A  (or  A  )  must  satisfy 
1  2  2  1 

the  following  relation  in  order  that  the  three-dimensional  theory  reduce 


to  uniaxial  Eq.  (110): 


•VT « -  ai> 


(120g) 
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Clearly,  there  is  a  single  damage  parameter  in  this  representation  which 

reflects  the  extent  of  microcracking;  it  is  weighted  Lebesque  norm 

||a,0  +  A  /?"  I  L  in  Eq.  (120f ) .  Moreover,  if  A  and  A  are  poeitive, 

1  2  2'  'nq  1  2 

pressure  (3  <  0)  is  seen  to  correctly  suppress  the  amount  of  crack  growth. 

It  should  also  be  recalled  that  M  reflects  the  effect  of  aging  and/or 
rehealing.  This  function  is  unity  in  the  absence  of  these  effects,  but 
otherwise  is  given  by 


M 


r 

r 

r 


(a  1 1  >  r 
m  1  r 


<°m  Il> 


(121) 


where  subscript  r  denotes  constant  quantities  corresponding  to  a  pre¬ 
selected  reference  state;  also,  recall  that  Dj  is  the  coefficient  in  creep 
compliance  Eq.  (93)  and  that  T  and  (o^  Ij)  are  fracture  properties. 

The  above  theory  was  developed  under  the  assumption  that  cracks  propa¬ 
gate  in  the  opening  mode.  However,  the  equations  governing  fracture  in 
shearing  modes  are  believed  to  be  anulogous  to  the  opening  mode  case  [3] 
and,  therefore,  we  suggest  that  Eq.  (120)  may  be  valid  regardless  of  the 
crack  mode.  However,  further  study  is  needed  to  confirm  this  point  or 
else  to  determine  if  the  more  general  theory,  Eq.  (47),  is  valid. 

It  is  important  to  recognize  that  the  three-dimensional  theory  based 
on  Eq.  (47)  along  with  Eq.  (120)  is  identical  to  linear  theory  except 
stresses  are  multiplied  by  a  scalar  factor,  l/aG.  Hence,  the  inverse  for¬ 
mulation  in  which  the  stress  tensor  is  expressed  as  a  functional  of  the 
strain  tensor  likewise  will  be  identical  to  linear  theory  except  is 

replaced  by  <J.,/a_.  This  latter  formulation  suffers  from  the  fact  that 
ij  G 

the  Lebesque  norm  is  expressed  in  terms  of  stress  invariants  rather  than 
strain  invariants;  consequently,  it  generally  will  not  be  possible  to 
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obtain  an  explicit  representation  for  stresses  in  terms  cr  strains  and 
strain  invariants. 

Let  us  now  consider  how  the  microcracking  model  cju  be  combined 
with  nonlinear  constitutive  equations  to  obtain  a  theory  which  is  poten¬ 
tially  applicable  to  a  wide  range  of  propellant  behavior.  First,  it 
should  be  recalled  that  Eq.  (110)  is  based  on  the  assumption  that  cracks 
are  (approximately)  arrested  after  rapid  growth.  As  the.  applied  stress 
is  increased  one  can  expect  these  previously  arrested  cracks  to  again 
propagate,  thereby  producing  significant  vacuole  dilatation  and  eventually 
gross  failure.  (It  should  be  added  that  significant  dilatation  does 
not  necessarily  result  from  the  first  growth  stage  when  the  binder  is 
relatively  soft;  e.g.,  when  a  crack  is  cut  into  a  sheet  of  rubber 
having  high  ultimate  elongation,  and  the  sheet  is  stretched  in  a  direc¬ 
tion  normal  to  the  original  crack  plane,  the  crack  may  open  so  much  that 
it  again  becomes  a  sharp  crack  with  its  plane  parallel  to  the  loading 
direction  [1].  The  net  change  in  dilatation  for  this  process  is  essen¬ 
tially  zero.)  If  it  can  be  assumed  that  the  subsequent  propagation 
obeys  the  previously  developed  model,  we  suggea,.  that  all  damage  may  be 
characterized  by  the  same  function  a^,  Kq.  :120e),  except  the  distribution 
function  G(g)  will  depend  on  strain  or  stri .  •.  invariants.  Moreover,  in 
the  high  stress  range  it  probably  will  be  necessary  to  account  for  non¬ 
linear  stress  invariant-dependence  in  Q  and  in  the  derivatives  of  F0 
appearing  in  Eq.  (47).  Similarly,  strain  formulation  Eq.  (30)  may  apply 
in  the  high  strain  range  after  the  stress  tensor  is  replaced  by  a^j/aG. 
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7.  Thermodynamic  Constitutive  Equations  with  Implicit  Pressure  Dependence 


The  strain  formulation,  Eqs.  (5)  and  (30),  and  special  cases,  e.g. 

Eqs.  (19)  and  (44),  are  based  on  the  assumption  that,  except  for  damage 
parameters,  the  nonlinearities  are  defined  by  strain  (rather  than  stress) 
invariants.  In  principle,  one  can  retain  the  relatively  simple  forms  of  the 
special  casas,  but  yet  account  for  large  superposed  pressures.  Referring  to 
the  thermodynamic  theory  in  [17],  it  is  easily  shown  that  the  Eqs.  (5)  and 
(30),  and  resulting  special  cases,  may  be  applied  to  deformations  superposed 
on  the  pressurized  state.  Namely,  one  interprets  the  stresses  and  strains, 
and  ,  as  values  which  are  added  to  the  values  existing  in  a  reference 
pressure  state.  All  material  properties  may  then  depend  on  this  pressure  as 
well  as  strain  invariants.  In  grain  analysis ,one  could  select  as  reference 
pressure  that  in  the  bore,  and  then  account  for  its  time-dependence  just  as 
one  does  for  a  spacewlse  uniform,  transient  temperature. 


G.  Generalization  for  Finite  Strains 


In  all  preceding  three-dimensional  equations  make  the  substitutions  [12], 

*  EU 

p  <122> 

po  ij 

0.,  -*■  -  T  J 

ij  P 

where  the  deviatoric  components  of  and  are  calculated  by  using  the 
same  definitions  as  for  small  strain  theory  and 

-  Lagrangian  (covariant)  strain  tensor 

■  (contravarlant)  stress  tensor  referred  to  convected 
coordinates  and  measured  per  unit  area  of  the  deformed 
body. 


p  ,  P  -  density  in  the  undeformed  and  deformed  states,  respec¬ 
tively. 


III.  EXPERIMENTAL  INVESTIGATIONS 


A.  Test  Equipment  end  Procedures 

Poker  chip  tests:  These  tests  were  made  on  ANB-3335-1  solid  propel¬ 
lant  samples  furnished  by  Aerojet  Solid  Propulsion  Company  (ASPC) .  As 
received,  the  samples  were  4.5  by  4.5  by  0.33  inches.  The  upper  and  lower 
surfaces  had  been  machined  to  attain  the  0.33  in.  thickness.  Therefore, 
before  attempting  to  bond  these  samples  to  the  plattens  used  for  loading, 
all  loose  oxidizer  particles  on  the  surface  of  the  samples  had  to  be  re¬ 
moved;  two  methods  were  tried.  The  first  method  was  to  spread  Duco  cement 
on  the  surface  of  the  propellant  sample  and  then  press  a  gauze  bandage  down 
firmly  onto  the  cement.  After  the  Duco  cement  dried  for  eight  hours  at 
room  temperature,  the  gauze  along  with  the  cement,  and  hopefully  the  loose 
oxidizer  particles,  were  gently  peeled  back  from  the  surface  of  the  pro¬ 
pellant  sample.  This  procedure  was  based  on  previous  experience  by  others, 
in  which  it  was  found  to  lead  to  very  good  bonds.  Two  initial  samples 
cleaned  by  the  above  technique  and  bonded  to  the  steel  plattens  turned  out 
to  fail  under  low  loads,  with  gross  tearing  occurring  near  the  bond  line 
where  the  material  was  extremely  soft  (probably  due  to  the  Acetone  in  Duco 
cement) . 

The  second  method  used  to  remove  the  loose  oxidizer  particles  was  to 
press  down  strips  of  three-’inch  wide  surgical  adhesive  tape  across  the 
surface  of  the  propellant  and  then  peal  back  gently.  All  samples  used  in 
the  mechanical  property  tests  were  prepared  in  this  way. 

Steel  plattens  0.25  in.  thick  and  4.00  inches  in  diameter  were  sand 
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blasted  by  a  medium  grit  In  order  to  prepare  them  for  bonding  to  the  sample. 
The  surfaces  were  then  cleaned  using  Acetone. 

Uralane  8615  epoxy  compound,  a  two-part  epoxy  mixed  In  approximately 
equal  parts  by  weight  and  obtained  from  Furane  Plastics,  Los  Angeles, 
California,  was  used  to  bond  the  propellant  samples  to  the  plattens.  This 
epoxy  is  a  black  material  which  has  a  pot  life  of  approximately  twenty 
minutes  once  the  two  components  are  mixed  and  a  complete  curing  time  of 
about  twenty-four  hours.  The  mixed  epoxy  was  placed  on  both  propellant 
surfaces  and  spread  as  thinly  as  possible  to  eliminate  bubbles.  Two  steel 
plattens  were  fixed  in  a  lathe  using  centered  pull  rods  which  had  been 
screwed  into  each  plate.  A  poker  chip  specimen  was  completed  and  made 
ready  for  curing  by  placing  the  propellant  coated  with  epoxy  between  the 
two  plattens  and  then  mechanically  forcing  the  plattens  against  the  samples. 
Since  the  samples  would  slide  out  if  left  alone,  several  pieces  of  tape 
were  used  to  hold  the  sample  centered  on  the  plattens.  The  bonded  sample 
was  allowed  to  cure  for  twenty-four  hours  in  this  position  at  room  tempera¬ 
ture. 

The  specimen  was  then  turned  in  a  lathe  at  a  very  slow  speed  and 
excess  propellant  was  trimmed  off  even  with  the  steel  plattens.  Completed 
poker  chip  specimens  were  stored  in  a  dry  refrigerator  maintained  at  74°F 
and  containing  a  desicant  to  maintain  low  humidity. 

Just  before  a  test,  a  teflon  coated  thermocouple  girth  wire  was 
placed  around  the  outside  of  the  specimen  midway  between  the  two  plattens, 
as  shown  in  Fig.  A.  One  end  of  the  girth  wire  was  stripped  and  two 
closely  spaced  dots  of  solder  were  placed  on  this  end  of  the  wire;  by 
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inserting  a  horseshoe-shaped  1/4  in.  Long  x  36  gage  wire  holder  in  between 
these  closely  spaced  solder  dots  and  into  the  propellant,  the  girth  wire 
was  secured  at  one  end.  This  wire  extended  from  the  tiedown  point  completely 
around  the  circumference  of  the  specimen  and  came  off  tangent  to  the  tie¬ 
down  point  to  then  connect  to  a  lever  arm.  Horseshoe-shaped  portions  of 
bare  36  gage  wire  were  stuck  into  the  propellant  over  the  girth  wire  at 
Intervals  of  approximately  0.4  inches,  to  hold  it  in  position  along  the 
mid-plane  between  the  two  plattens. 

An  aluminum  plate  approximately  1/4  in.  thick  by  6.5  in.  diameter 
was  mounted  on  top  of  the  upper  platten  and  fastened  securely  with  a  nut 
on  the  pull  rod.  This  upper  plate  was  designed  to  accommodate  three  ver¬ 
tical  micrometers  at  120  degree  spacings.  Another  aluminum  plate  of  the 
same  size  and  an  S-shaped  bracket  were  mounted  securely  on  the  lower  side 
of  the  lower  platten  by  tightening  the  nut  on  the  lower  pull  rod.  This 
lower  aluminum  plate  was  designed  to  accommodate  three  linear  variable 
differential  transformers  (LVDT's)  aligned  with  shafts  from  the  three 
vertical  micrometers  on  the  upper  plate.  These  vertical  micrometers  had  a 
total  travel  of  0.05  in.  which  allowed  us  to  make  very  fine  in-place  cali¬ 
brations. 

For  tension  tests  the  upper  pull  rod  was  attached  to  a  universal 
Joint  which  was  in  turn  attached  to  a  single  lever  arm  loading  frame, 

(see  Fig.  5).  The  lower  pull  rod  was  attached  to  a  universal  joint  which 
was  in  turn  attached  to  a  slip  joint  which  engages  only  when  the  lever  arm 
is  loaded.  An  LVDT  was  attached  to  the  S-shaped  bracket  that  extends  from 
beneath  the  lower  platten  to  be  in  line  with  the  girth  wire.  The  other  end 
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Figure  5.  Poker-Chip  Assembly  used  for 
Tensile  tests. 
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of  the  LVDT  was  attached  to  a  Chavetz  micrometer  which  had  a  total  travel 
of  1.0  in.  A  modified  4-4  1/2  in.  hose  clamp  was  attached  to  the  lower 
platten  to  support  a  small  L-ahaped  lever  arm.  The  vertical  arm  of  the 
lever  was  connected  to  the  girth  wire,  while  the  horizontal  arm  was  con¬ 
nected  to  a  vertical  wire  extending  to  a  slug  within  the  LVDT  mounted  on 
the  lower  S-6haped  bracket.  The  horizontal  girth  motion  was  thus  trans¬ 
ferred  into  a  vertical  motion  and  magnified  by  approximately  four.  With 
the  Chavetz  micrometer  it  was  possible  to  attain  an  inplace  calibration 
for  the  girth  measurement.  To  insure  that  there  was  no  sticking  associated 
with  the  girth  wire  or  LVDT's,  a  pulsing  vibrator  was  attached  to  the  frame 
of  the  test  mechanism. 

The  three  vertical  LVDT's  spaced  120  degrees  apart  around  the  sample 
were  used  to  measure  plate  separation  and  to  determine  if  there  was  bend¬ 
ing  in  the  sample.  By  averaging  the  three  readings  we  obtained  an  overall 
average  strain  which,  together  with  the  girth  measurement,  enabled  the 
determination  of  dilatation.  Signals  from  all  four  LVDT's  were  fed  into 
an  amplification  system  that  allowed  us  to  obtain  high  gain  settings  and 
to  maintain  a  linear  output  signal  as  a  function  of  displacement  for  a 
given  range. 

For  testing,  the  system  was  loaded  to  provide  10,  20,  35,  and  50  psi 
average  axial  stress  to  the  poker  chip  specimen.  Testing  was  performed  at 
a  temperature  of  74°F  and  at  a  relative  humidity  of  between  34  and  40  per¬ 
cent.  All  tests  were  creep  and  recovery  tests  with  cycles  of  five  minutes 
creep,  ten  minutes  recovery,  and  15  minutes  wait  between  cycles.  There 
were  one  or  two  one-minute  initial  loadings  at  each  load  range  before  the 
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actual  teat  to  Insure  that  the  data  readings  would  remain  on  scale  and  that 
the  readings  were  still  linear.  Five  to  seven  creep  and  recovery  cycles 
were  run  for  each  load  range,  and  only  one  load  range  was  run  per  day.  Com¬ 
pression  tests  were  also  run  on  poker  chip  specimens.  The  only  difference 
in  the  compression  test  assembly  was  that  Instead  of  having  the  upper  and 
lower  universal  joints  attached  to  the  upper  and  lower  pull  rods,  respectively, 
heavy  steel  frames  replaced  the  pull  rods;  one  frame  was  used  to  attach  the 
lower  platten  to  the  upper  universal  joint,  and  another  frame  connected  the 
upper  platten  to  the  lower  universal  joint. 

Strip-biaxial  teats:  Tests  were  performed  on  ANB-3066  and  ANB-3335-1 
solid  propellant  obtained  from  Aerojet  Solid  Propulsion  Company.  The 
biaxial  dilatometer  used  for  these  tests  was  obtained  from  Hill  Air  Force 
Base  on  a  loan  basis  and  as  received  was  in  unworkable  condition  as  originally 
manufactured  by  the  CETEC  Corporation.  The  unit  was  completely  dissassembled , 
and  after  extensive  modification  was  put  into  working  condition.  The  dila¬ 
tometer,  shown  in  Fig.  6,  was  designed  to  measure  dilatation  by  two  differ¬ 
ent  methods.  The  first  method  relied  upon  direct  LVDT  measurements  of  decrease 
in  cross-sectional  thickness.  The  LVDT  system  consisted  of  spring  loaded 
probes  that  slide  on  the  surface  of  the  specimen's  sides.  Unfortunately 
the  probes  were  not  only  insensitive  to  the  small  thickness  changes  but  also 
damaged  the  surface  of  the  propellant  at  the  point  of  contact.  As  a  result 
the  LVDT  probes  were  removed  entirely,  and  we  relied  upon  the  second  method 
of  dilatation  measurement.  This  method  was  based  on  state-of-the-art  gas 
dilatometric  techniques  such  as  those  employed  by  Farris  [11]. 

The  biaxial  dilatometer  consists  of  two  chambers,  a  specimen  chamber 
and  reference  chamber,  which  are  completely  separate  from  each  other.  The 


Figure  6.  Biaxial  Di.latometer  Mounted  in 
instron  Test  Machine. 
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dilatometer  was  configured  so  thaC  both  chambers  could  be  pressurized  up 
to  1000  psi  and  equalized  to  the  same  pressure.  Once  all  the  leaks 
around  the  dilatometer,  dumbing  connections,  door,  shaft,  and  calibration 
micrometer  were  sealed  at  high  pressure,  the  system  was  functional.  The 
dilatometer  was  mounted  in  an  Instron,  and  an  extensometer  was  mounted  to 
the  shaft  of  the  dilatometer  to  measure  sample  extension.  Unlike  other 
dilatometers  which  have  a  compensating  shaft  that  extends  from  the  reference 
chamber,  the  dilatometer  used  in  these  tests  was  not  equipped  with  this 
reference  chamber  shaft;  specimen  chamber  shaft  movement  was  compensated 
electronically  by  coupling  the  shaft  extension  with  the  pressure  change 
due  to  shaft  movement  alone. 

Before  each  test,  the  unit  was  calibrated  for  volume  change  by 
inserting  a  precision  micrometer  shaft  into  the  specimen  chamber  in 
increments  to  simulate  a  volume  change  of  .00707  in  (for  the  samples  used 
in  our  tests  this  corresponds  to  a  volume  change  of  about  0.5  percent). 

It  was  possible  to  resolve  volume  changes  on  the  order  of  .01  to  .02 
percent  over  a  relatively  long  period,  provided  the  temperature  remained 
constant.  To  insure  constant  test  conditions,  the  tests  were  run  in  large 
walk-in  chambers  which  were  held  to  a  constant  temperature  and  humidity. 

The  dilatometer  and  associated  plumbing  were  wrapped  thoroughly  with  insula¬ 
tion  material  (see  Fig.  6)  to  prevent  changes  due  to  the  body-heat  of  people 
in  the  room  who  are  running  the  tests. 

The  propellant  samples  furnished  by  Aerojet  were  saw-cut  into  speci¬ 
mens  with  dimensions  of  7.0  by  1.63  by  0.38  inches.  Using  a  circular  cutter, 
each  end  of  each  specimen  was  trimmed  to  remove  the  material  included  within 
a  circular  arc  with  an  approximately  one-inch  radius  connecting  the  two  corners 


at  each  end  of  the  specimen,  as  shown  in  Fig.  7.  The  samples  were  then 
bonded  to  aluminum  plates  using  the  8615  epoxy  described  in  the  poker- 
chip  section  above.  The  aluminum  plates  were  prepared  by  first  sand 
blasting  the  surface  to  be  bonded  and  then  placing  masking  tape  around  the 
plates  so  that  the  tap«.  extended  approximately  0.25  inches  above  the  sur¬ 
face  of  the  plates.  The  epoxy  was  then  poured  into  the  cups  made  by  the 
masking  tape  to  about  half  full.  One  7.0  by  0.38  inch  edge  of  the  sample 
was  pushed  into  the  epoxy  and  against  the  aluminum  plate.  Care  was  taken 
to  create  a  fillet  around  the  sample  as  it  emerged  from  the  epoxy  in  order 
to  provide  for  gripping  and  a  gradual  load  transfer  from  the  grips.  A 
spacer  was  placed  on  each  side  of  the  samples  to  hold  the  propellant  at  the 
center  of  the  aluminum  plate  while  the  epoxy  cured  for  twenty-four  hours. 
This  process  was  then  repeated  for  the  opposite  edge  of  the  sample.  After 
curing,  the  bonded  samples  were  stored  for  at  least  twenty-four  hours  be¬ 
fore  testing  in  constant  temperature  and  humidity  walk-in  chambers. 

A  bonded  sample  was  fastened  to  loading  rails  and  loaded  info  the 
dilatometer  specimen  chamber  as  shown  in  Fig.  8.  The  door  was  then  put 
in  place  and  bolted  down. 

If  the  test  was  to  be  run  under  pressure,  the  valve  connecting  the 
reference  and  specimen  chambers  was  opened  and  the  system  was  pressurized 
to  the  desired  pressure.  Thu  valve  was  then  closed  to  isolate  the  cham¬ 
bers  again  while  the  system  stabilized  for  one  to  two  hours.  Just  before 
running  the  actual  test,  the  two  chambers  were  equalized  by  opening  the  con¬ 
necting  valve  and  then  closed  to  isolate  them  again.  Also,  as  described 
above,  the  unit  was  calibrated  for  volume  change.  In  addition,  the  internal 
load  cell  located  within  the  dilatometer  was  calibrated. 
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Figure  7.  Unfailed  and  Failed  Strip-biaxial 
Specimen  End-Bonded  to  Aluminum. 


Several  types  of  tests  were  run  Including  constant  strain  rate  to 
failure,  ramp-relaxation,  and  cyclic  loading  and  unloading  before  load¬ 
ing  to  failure.  Sample  displacement,  volume  change,  and  ,ioad  were  re¬ 
corded  on  a  light  sensitive  oscillograph  recorder. 

Uniaxial  tests:  These  tests  were  run  on  AND-3066  and  ANB-3335-1 
Aerojet  propellant  using  the  same  biaxial  dilatometer  used  for  the  strip- 
biaxial  tests.  The  uniaxial  samples  were  made  from  the  7.0  by  1.63  by 
0.38  inch  strip-biaxial  specimens.  Instead  of  cutting  circular  arcs  of 
material  away  from  the  ends  as  was  done  for  the  biaxial  samples,  the 
rectangular  samples  were  cut  completely  through  to  within  0.2  inches  of 

each  edge  at  intervals  of  0.37  inches  along  the  length  of  the  sample.  The 
samples  were  then  bonded  to  the  aluminum  plates  as  described  for  biaxial 
specimens.  This  procedure  provides  nineteen  uniaxial  samples  to  be  tested 
at  one  time,  which  are  shown  in  Fig.  9. 

B.  Experimental  Data  and  Comparison  with  Analysis 

Poker-chip  tests:  Figures  10-13  show  compliances  and  dilatation  ob¬ 
tained  from  tension  tests  performed  on  samples  Nos.  1  and  2,  while  Figs. 
14-17  show  these  quantities  for  compression  of  sample  No.  2.  The  creep 
and  recovery  compliances  are  defined  as  follows: 

D  s  £  (123a) 

pc  a 

D  (r)  =  —  (123b) 

pc  a  v  ' 

where  o  is  the  total  axial  load  divided  by  platen  area,  e  is  the  strain 

(based  on  platen  separation)  during  the  creep  portion  of  one  loading  cycle. 


A-63 


I  ill  rrq - 1  i 1 — |  i  vr-q - r— i  rn 


B- 


_  _  b a - “ — 1 fl~' h 


Figure  10.  Poker  Chip  Creep  Compliance  in  Tension. 
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Figure  13.  Poker  Chip  Dilatation  in  Tension-Recovery 
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Figure  14.  Poker  Chip  Creep  Compliance  in  Compression. 
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and  er  is  the  recovery  strain  (which  is  the  strain  following  load  removal). 
The  cycle  numbers  shown  on  the  graphs  represent  a  loading  and  unloading 
sequence,  with  cycle  No.  1  defining  the  first  sequence  after  the  specimen 
had  been  loaded  once  or  twice  for  calibration  purposes.  Therefore  some 
damage  had  alread  ie.n  incurred,  even  before  the  data  labeled  cycle  No.  1 
was  obtained.  N- jowever,  that  in  most  cases  for  each  stress  level  no 
further  significant  increase  in  damage  (as  measured  by  changes  in  compliance 
and  dilatation; appears  to  have  occurred  aa  the  sample  was  cycled. 

The  dilatation  ahovm  in  these  figures  was  found  using  the  same  pro¬ 
cedure  described  in  [  4]  for  earlier  work  at  TAMU;  specifically,  the 
following  equation  was  used: 


V  ,  -  V 


where 


=  volume  displaced  by  moving  platen 

V  =  volume  swept  out  by  revolving  the  parabolic  area 
p  at  the  periphery  around  the  loading  axis. 

As  an  aid  in  interpreting  the  compliance  and  dilatation  data,  we  show 
in  Fig.  18  a  graph  that  enables  various  property  ratios  tc  be  deduced  from 
the  experimental  data.  This  graph  was  obtained  from  a  linear  quasi-elastic 
analysis  of  a  compressible  poker-chip  using  the  equations  in  [  4 ) ;  a  con¬ 
densed  version  of  this  graph  also  appeared  in  [  4 ) .  The  following  definitions 


apply: 


Epc  =  Dpc  -  effective  relaxation  modulus  of  poker  chip 


"  creep  compliance  under  uniaxial  loading 
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B  ■  bulk  creep  compliance 


K  »  bulk  relaxation  modulus  (=  B 

/  -1s 

£  ■  relaxation  modulus  under  uniaxial  loading  (-  D  ) 

u  u 

In  order  to  illustrate  the  use  of  this  graph  let  us  consider  some  of 

the  data  in  Figs.  10-17.  The  ratio  AV/Vqe  is  reasonably  constant  in  Fig.  12 

and,  for  sample  No.  2,  has  the  approximate  value  of  0.7.  Reading  across 

and  down  from  the  AV/V  e  curve  in  Fig.  18  we  find  D  /D  a  E  /E  “3.5. 

o  °  u  pc  pc  u 

Reading  up  to  the  B/D  curve  and  over  to  the  ordinate  we  find 
pc 

B/D  “  E  /K  “  0.9. 
pc  pc 

From  the  compliance  curve  in  Fig.  10  one  can  obtain  D  c  at  any  time, 
and  therefore  calculate  the  remaining  properties  in  these  ratios.  For 
example,  at  t  ■  0.1  min,  D  ^  ~  3000  psi,  and  therefore  K  “  3300  psi.  This 
extremely  low  value  predicted  for  bulk  modulus  (which  is  typical  of  the 
values  obtained  from  the  other  data)  is  probably  in  error  because  of  non¬ 
linear  effects.  However,  it  is  believed  that  its  order  of  magnitude  is  not 
unreasonable  since  the  specimen  probably  has  a  significant  void  fraction  as 
a  result  of  prior  loadings. 

An  observation  of  interest  is  that  if  Av/V^e  is  independent  of  time, 
as  it  is  in  Fig.  12  and  in  Fig.  16  except  at  short  times,  linear  theory 
can  be  used  to  show  that  Poisson's  ratio  is  constant.  Now,  a  composite 
material  consisting  of  a  relatively  incompressible  viscoelastic  binder  to¬ 
gether  with  voids  and  rigid  particles  can  be  shown  to  have  a  constant  Poisson's 
ratio;  this  result  is  easily  established  by  means  of  the  linearity  assumption 
plus  dimensional  analysis.  This  latter  observation  is,  of  course,  consistent 
with  the  hypothesis  that  the  poker  chip  specimens,  as  tested,  contained  a 
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large  initial  void  fraction  due  to  previously  induced  damage. 

A  partial  check  on  the  microcrack  theory  in  Section  II-E  1b  pro¬ 
vided  by  a  comparison  of  creep  nnd  recovery  data.  First  of  all,  conaider 
the  following  Isothermal  constitutive  equation  relating  average  axial 
stress  applied  to  the  poker-chip,  and  the  axial  strain,  c: 

e  -  /  D  (t  -  f)  dt'  (125) 


where  P  -  P(o)  is  a  nonlinear  function  of  stress.  Also  is  the  linear 
viscoelastic  poker-chip  creep  compliance.  The  creep  and  recovery  com¬ 
pliances  are  predicted  to  be,  respectively, 


D  (t)  -  -  D  (t) 
pc  a  p 

D  (r)  “  £  [D  (t)  -  D  (t  -  t  )] 
pc  aip'7  pv  c'J 


(126a) 


(126b) 


where  t£  is  the  time  at  which  the  load  is  removed.  Combining  Eqs.  (126a) 
and  (126b)  we  find 

D  ^  ■  D  (t)  -  D  (t  -  t  )  (127) 

pc  pc  pc'  c' 

Since  Eq.  (125)  with  P  -  a  is  an  exact  equation  for  linear  materials,  we 
conclude  that  the  recovery  compliance  for  this  special  nonlinear  theory  is 
related  to  creep  compliance  in  the  same  way  as  in  linear  theory.  Equation 
(127)  has  been  used  to  construct  the  theoretical  recovery  curves  shown  in 
Figs.  11  and  15.  Considering  the  scatter  in  experimental  data,  predictions 
for  the  tension  test  are  quite  good.  However,  the  theory  is  seen  to  under¬ 
predict  recovery  compliance  in  compression;  the  difference  is  believed  to  be 
due,  at  least  in  part,  to  rubbing  between  flaw  faces  which  would  Impede 
recovery  from  a  compressed  state.  In  both  tension  and  compression 
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the  creep  compliance  varies  appreciably  with  stress  level. 

Now,  refer  to  the  theory  with  microcracking.  Rigorously,  a  nonlinear 
stress  analysis  of  the  poker-chip  is  required  to  interpret  the  data  be¬ 
cause  the  stresses  are  not  spacewise  uniform.  However,  we  will  argue  that 
Eq.  (110)  is  approximately  valid  for  the  tension  teats  after  replacing  D 
by  0^,  if  the  underlying  three-dimensional  theory,  Eqs.  (47)  and  (120), 
are  themselves  correct. 

Considering  linear  theory  first,  if  Poisson's  ratio,  v,  is  constant 
and  the  material  is  linear,  the  stresses  throughout  the  poker-chip  can  be 
shown  to  be  timewise  constant  during  creep  and  zero  during  recovery.  In 
tension  tests,  the  thickness-averaged  normal  stresses  are  approximately 
equal  and  distributed  parabolically  across  the  radius  (with  the  maximum 
occurring  at  the  center),  while  the  thickness-averaged  shear  stress  in¬ 
creases  linearly  with  radius  from  the  center  [18].  With  this  observation 
in  mind,  and  assuming  the  stresses  in  a  nonlinear  poker-chip  are  dis¬ 
tributed  in  the  same  general  way,  we  suggest  that  the  Lebesque  norm  of 
stress  invariants  in  Eq.  (120f)  may  be  reasonably  uniform.  On  physical 
grounds  it  can  be  argued  that  A^  and  A^  are  positive;  therefore  the 
variation  in  0  will  tend  to  cancel  that  of  since  the  former  invariant 
is  large  when  the  latter  one  is  small  and  vice-versa.  (This  cancellation 
of  effects  will  not  occur  with  compression  since  both  invariants  increase 
algebraically  with  radius.)  If  indeed  g',  and  therefore  a^,  are  essen¬ 
tially  independent  of  radius,  the  relation  between  average  stress  and 
axial  strain  will  be  essentially  the  same  as  for  linear  theory  except 
a/a^  will  replace  stress;  this  can  be  shown  by  referring  to  the  thickness- 
averaged  field  equations  for  the  poker-chip  [19). 
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Hence,  by  making  the  substitution  P  -*■  o/aG  in  Eq.  (125),  the 
approximate  equation  with  microcracking  is  obtained.  That  this  equation 
successfully  predicted  recovery  from  creep  would  seem  to  imply  aG  is 
essentially  independent  of  time,  although  it  depends  on  stress  in  view 
of  the  observed  stress-dependence  of  creep  compliance.  For  the  power- 
law  distribution  function  Eq.  (Ill),  we  find  during  creep  (0  <  t  <  t  ) : 

r-1 

7"  “  1  +  C  (a)r_1(t)q  (128) 

G 

where  C,  is  a  constant;  during  recovery  a-  ^  is  given  by  £q.  (128)  with 

1  Ur 

t  -  t^.  For  the  data  in  Fig.  10,  n  =  1/6,  which  implies  q  «  14.  More¬ 
over,  we  find  the  stress-dependence  of  the  creep  compliance  in  Fig.  10 
can  be  predicted  by  neglecting  the  unity  term  ir  Eq,  (128)  and  setting 
r  “  1.8.  Thus,  the  time-dependence  of  aG  is  relatively  weak  compared 
to  that  of  the  creep  compliance  itself.  It  is  of  interest  to  note  in 
passing  that  the  value  r  -  2  is  predicted  from  an  approximate  distribution 
function  for  strain  concentration  factors  derived  by  Farris  [19 ]  for 
spherical  particles. 

Strip-biaxial  and  uniaxial  tests:  A  large  number  of  constant  strain 
rate  tests  were  conducted  at  room  temperature  and  at  approximately  -20°F 
in  order  to  assess  the  accuracy  of  simplified  constitutive  equations  in 
Sections  II-A  and  II-D. 

Although  most  of  the  tests  were  conducted  on  ANB-3335-1  propellant, 
a  study  of  ANB  3066  at  T  *  -20°F  was  made  early  in  the  program;  typical 
results  for  stress  and  dilatation  are  shown  in  Fig.  19.  The  prediction 
of  dilatation  was  made  by  first  evaluating  the  constants  ag  and  bg  in  the 
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LATATION  AV/V  (*) 


MSS  theory  and  a  and  b  in  the  OSS  theory  from  the  uniaxial  data.  Equa¬ 
tions  (86b)  and  (88b)  were  then  used  to  construct  the  dilatation  for  the 
biaxial  teat.  Obviously,  the  octahedral  strain  theory  is  moat  accurate 
for  these  data.  Also,  referring  to  Eq.  (32),  we  find  if  0  “  2  the  stress 
can  be  predicted  reasonably  well. 

Similar  results  were  found  for  ANB  3335-1,  as  shown  in  Fig.  20.  Again 
the  OSS  theory  is  seen  to  be  better  than  the  MSS  theory.  Also,  for  this 
case  0  =<  2.  Typical  room  temperature  data  on  this  propellant  are  given  in 
Fig.  21,  for  which  0  “  2.5.  The  biaxial  dilatation  data  now  fall  between 
the  OSS  and  MSS  theories.  A  study  of  the  data  obtained  on  several  other 
specimens  . evealed  that  the  OSS  prediction  was  generally  best,  although 
some  of  the  results  were  similar  to  that  shown  in  Fig.  21.  In  most 
of  these  latter  cases  some  of  the  uniaxial  specimens  cut  from  a  single 
strip  appeared  to  fail  prematurely,  which  would  account  for  the  OSS  pre¬ 
diction  being  low  compared  to  the  biaxial  dilatation  data. 

In  order  to  check  the  hypothesis  that  I3  is  absent  from  the  dilata- 

tional  constitutive  Eq.  (10)  (or  Eq.  (61)),  it  1b  not,  of  course,  necessary 

to  use  a  piecewise  linear  representation  of  dilatation.  If  the  term  0/3Ke 
is  negligible  in  these  equations,  then  a  plot  of  the  actual  uniaxial  and 
biaxial  data  on  a  graph  of  1^  vs.  will  form  a  single  curve  if  I3  does 
not  enter  the  constitutive  equation.  Figures  22-25  show  such  plots  for 
0  pel  and  200  psi.  Although  there  is  some  spread,  the  theory  seems  to  be 
quite  adequate  for  correlating  the  dilatation  in  the  biaxial  and  uniaxial  modes; 
most  of  the  spread  is  believed  due  to  premature  failure  of  some  uniaxial  samples. 

It  should  be  added  that  in  view  of  the  relatively  large  strains  involved, 

we  have  used  invariants  based  on  Lagrangian  strains  in  Figs.  22-25;  see 


STRAIN  (%) 

Figure  20.  Uniaxial  and  Strip-Biaxial  Data  on  ANB  3335-1  Pr< 


Figure  21.  Uniaxial  and  Strip-Biaxial  Data  on  ANB  3335-1  Propellant. 
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T  -  -16°F 
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Figure  22.  Relation  Between  First  and  Second  Strain  Invariants  for 

Data  in  Fig.  20. 
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Figure  23.  Relation  Between  First  and  Second  Strain  Invariants 
for  ANB  3335-1  Propellant 
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STRAIN  RATE 
(ml n-1 ) 

Unlaxi.;! 

0.5 

Biaxial 

0.5 

Uniaxial 

0.5 

Biaxial 

0.5 

Figure  24.  Relation  Between  First  and  Second  Strain  Invariants 
for  ANB  3335-1  Propellant 
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Figure  25.  Relation  Between  First  and  Second  Strain 
Invariants  for  ANB  3’35-l  Propellant. 


Section  II-G,  la  particular,  with  primes  denoting  daviatoric  strains, 


IT  H  Btl  -  3aAT  (129e) 

V  V  <129b> 

Note  that  the  invariant  defined  by  Eq.  (129a)  ia  not  equal  to  the  dilata¬ 
tion,  except  for  email  atraina. 

Finally,  it  ie  of  interest  to  compare  the  value  of  the  modulus  ratio 
in  Eq.  (23)  found  from  the  vriaxlal  and  biaxial  stress-strain  data  with 
that  obtained  from  the  poker  chip  tests.  For  the  former  tests  at  74°F  we 
have  found  that  the  relation  S/a  »  5  holds.  Thus, 


«  5/S  “  12  (130) 

e 

Referring  to  Fig.  18,  it  is  found  that  this  result  predicts  the  dilatation- 
strain  ratio 

AV 

-  0.69  (131) 

where  the  standard  linear  elastic  relation 


3K  E 

G  -  e~u 
e  9K  -  E 


(132) 


has  been  used;  also,  corresponding  to  the  notation  in  Fig.  18,  K  =  K.  The 

6 

ratio  in  Eq.  (131)  is  of  the  correct  order  of  magnitude,  based  on  the  poker 


IV.  CONCLUSIONS 


itephaals  in  this  study  has  been  on  (i)  establishing  a  simple  consti¬ 
tutive  theory  for  propellant  at  large  values  of  vacuole  dilatation ,  and 
(li)  the  use  of  viscoelastic  fracture  mechanics  as  a  means  of  relating 
overall  mechanical  response  to  fundamental  parameters  that  define  the 
extent  of  microstructural  damage. 

Under  itam  (1)  we  hava  found  the  dilatation  of  strip-biaxial  speci¬ 
mens  and  uniaxial  specimens  to  be  essentially  ooual  when  referred  to  equal 
values  of  the  octahedral  ihear  strain.  An  implication  of  this  result  is 
that  the  third  devl&torit  strain  invariant  probably  can  be  omitted  from 
the  nonlinear  constitutive  theory  without  introducing  appreciable  error. 

It  is  shown  that  dil  .tat ion  measured  on  thin  simple  shear  specimens  would 
provide  another  check  of  this  conclusion;  however  shear  data  were  not  ob¬ 
tained  in  this  study. 

In  the  area  (il)  we  developed  a  constitutive  theory  which  is 
Identical  to  linear  viscoelasticity  except  the  stress  tensor  is  multi¬ 
plied  by  a  scalar  function  that  reflects  microstructural  damage.  This 
function  depends  on  a  weighted  Lebesque  norm  involving  stress  invariants 
(rather  than  strain  invariants)  and  is  expressed  in  terms  of  temperature- 
reduced  time;  the  weighting  function  in  the  norm  reflects  one  effect  of 
aging  and  rehealing,  if  any.  The  fracture  mechanics  theory  was  developed 
during  the  latter  portion  of  the  program  and  therefore  only  a  limited 
check  of  the  theory  is  given. 
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APPENDIX  B 


PREPROCESSOR  CODE  NLOOI 


1 .0  INTRODUCTION 

NLOOI  Is  a  preprocessor  for  the  nonlinear  viscoelastic  characteri¬ 
zation  code  NL002.  It  takes  the  raw  test  data,  consisting  of  a  temperature, 
strain  rate,  time  step  Increment,  stress  and  dilatation,  and  then  calculate 
the  following  data: 

1.  Principal  strains,  shear  strain  and  true  stress 

2.  Strain  rates 

3.  Strain  Invariants 

4.  Dilatation,  corrected  for  pressure  and  temperature 

5.  Octahedral  strain 

6.  The  10th,  20th,  30th,  40th,  50th,  60th  and  infinite  Lebesque  norms 

All  these  data  are  generated  for  the  same  time  points  as  the  original  data. 

The  generated  and  raw  test  data  are  then  assigned  a  user  prescribed  Identifi¬ 
cation  and  written  out  onto  a  tape  as  a  single  element.  This  process  is  then 
repeated  once  for  each  test,  with  each  test  beina  limited  to  100  input  data 
points.  There  Is,  however,  no  limit  to  the  number  of  tests  which  may  be  pro¬ 
cessed  in  a  single  run. 

Reproduced  from 
best  available  copy. 
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2.0 


PROGRAM  INSTRUCTIONS 


The  purpose  of  this  section  Is  to  list  the  basic  variables  of  the 
NL001  code  and  provide  a  set  of  input  Instructions  for  running  the  program. 

2.1  BASIC  VARIABLES 


Following  Is  a  list  of  the  basic  program  variables  of  NL001 . 
Each  variable  Is  described  briefly  and  labeled  as  either  an  input  or  calcu 
lated  output  variable.  All  Input  Is  also  saved  as  output.  Any  variable 
which  Is  not  either  input  or  output  Is  merely  a  working  variable  used 
Internally,  usually  for  Intermediate  calculations. 


Al,  A2 
BETA 
BULK 
Cl,  C2 


Working  variables 

Volumetric  expansion  coefficient  -  Input 
Bulk  modulus  -  input 

Arbitrarily  small  constants  -  fixed  In  program  by 
user.  Presently  1 .OE-6. 


OIL 

DR1,  DR2,  DR12 
DT 

EL'CT 

EPS 

E11 »  E22»  E33 

El  2 

INV1 

INV2 


Dilatation  an  input  -  corrected  for  temperature  and 
pressure  for  output 

Normal  and  shear  calculated  strain  rates  -  output 

Time  step  in  minutes  -  input 

Octahedral  strain  -  output 

Epsilon  tester  working  variable 

Normal  strains  -  output 

Shear  strain  -  output 

First  strain  invariant  =  +  E22  +  E33  -  output 

Second  strain  invariant  =  E^  E22  +  E22  E33  +  E33  En 
-  output 
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INV3 

KODE 


KTEMP 

LAST 


MATID 

NAME 

NDP 

NJ9RMF 

NORM1 

PRES 

RATE 

SP 

STRESS 

STRUE 

T 

TEMP 


Third  strain  Invariant  ■  E,2  E33  -  output 

Input  control  variable:  “  1  If  a  uniaxial  test, 

■  2  If  a  biaxial  test, 

■  3  If  a  shear  test 

Input  control  variable:  ■  1  If  a  constant  temperature  test 

■  2  If  a  variable  temperature  test 

Input  delimiter  flagging  last  Input  data  point. 

Set  blank  or  zero  if  not  last  data  point  -  nonzero  If 
last  data  point. 

Material  Identification  -  input 

Test  Identification  -  Input 

Number  of  data  points  of  a  specific  test 

Infinite  Lebesque  norm  -  output 

1th  Lebesque  norm  -  output 

Test  pressure  -  Input 

Strain  rate,  observed  -  input 

Working  variable 

Observed  stress-input 

True  stress  -  output 

Time  In  minutes  -  output 

Test  temperature  -  input 


2.2  INPUT  INSTRUCTIONS 

The  following  series  of  cards  are  read  once  for  each  test  processed 
through  NL001 .  The  last  card  of  the  input  deck  is  always  a  blank,  however  there 
is  no  blank  card  between  tests.  The  required  input  variables  fo**  each  card 
appear  below  with  their  formats  shown  in  parentheses. 
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Card  1  (A6,  E14.4,  215,  5a4,  2E15.6)  NAME,  PRES,  KTEMP ,  K0DE,  MATID,  BULK,  BETA 
Card(s)  2  (5E10.3,  15)  TEMP{J),  RATE(J),  DT(J),  STRESS(J),  OIL(J),  LAST 

Here  J  Indicates  the  Jth  data  point  of  this  test.  Card  2  is  repeated 
once  for  each  data  point. 

Cards  1  and  2  are  repeated  once  for  each  test  processed. 

Card  3  -  A  blank  card  indicating  that  last  test  has  been  processed. 

2.3  SAMPLE  INPUT  SHEET 

Shown  below  is  a  sample  Input  sheet  for  NL001 .  It  contains  the 
Input  data  for  two  uniaxial  constant  temperature  tests.  The  first,  U00100 
is  at  ambient  pressure  while  U10105  Is  run  at  100  pslg. 

2.4  PROGRAM  LISTING 

The  listing  of  NL001  appears  below.  Note  that  It  is  written 
specifically  for  the  UNIVAC  1108  computer.  Use  of  this  program  on  any 
other  computer  would  require  substitution  of  the  appropriate  statements 
for  the  *'AKELT  and  SQUEEZ  subroutines. 
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APPENDIX  C 


POST  PROCESSOR  CODE 


1 .0  INTRODUCTION 

POST  Is  a  postprocessor  for  the  NL001  code.  It  takes  the  output  of 
NL001  off  the  NL001  output  tape,  and  prints  selected  output  In  a  report 
format.  There  Is  a  minimum  of  two  pages  of  output  for  each  test  reported. 

The  first  two  columns  of  each  page  are  the  data  point  numbers  and  time  In 
seconds.  Page  1  then  contains  the  three  strain  Invariants,  the  dilatation 
aV/Vq,  corrected  for  any  temperature  and  pressure  effects,  and  the  octahedral 
strain.  Page  2  displays  the  temperature,  calculated  strains,  and  the  true 
stress.  In  addition  to  the  above,  each  page  has  a  heading  giving  material 
type,  bulk  modulus  and  volumetric  expansion  coefficient,  the  test  type, 
Initial  strain  rate,  pressure,  and  whether  It  Is  a  constant  or  variable 
temperature  test. 

The  only  card  input  to  POST  is  the  last  card,  a  blank,  signifying 
the  end  of  the  data.  All  other  Input  Is  read  directly  from  the  NL001  output 
data  tape. 

2.0  PROGRAM  INSTRUCTIONS 

The  purpose  of  this  section  is  to  list  the  basic  variables  of  the  POST 
code  and  provide  a  set  of  Instructions  for  running  the  program. 

2.1  BASIC  VARIABLES 

Following  is  a  list  of  the  basic  program  variables  of  POST.  The 
names  and  meaning  of  variables  Is  compatible  with  the  NL001  code. 


BETA 

V''umetr1c  expansion  coefficient 

BULK 

Bulk  modulus 

DUM 

Dumn\y  variable 

ED0TO 

Initial  strain  rate 

ITEMP 

Alphanumeric  constant/variable  temperature  header 

ITYPE 

Alphanumeric  unlaxlal/biaxial/shear  test  type  header 

K0DE 

Test  type  control  variable 

KTEMP 

Temperature  control  variable 

LAST 

Last  data  point  delimiter 

MATID 

Material  Identification 

NAME 

Test  element  identification 

NDP 

Number  of  data  points  in  a  specific  test 

PRES 

Test  pressure,  psig 

X 

An  array  used  to  store  the  various  output  quantities 

2.2  INPUT 

INSTRUCTIONS 

There  is  no  card  input  for  P0ST.  All  input  Is  read  from  the 
NL001  data  output  tape,  element  by  element.  The  user  merely  Inserts  an 
"ADD"  card  tor  each  element  of  the  tape  he  wishes  processed  by  POST.  After 
the  last  "ADD"  card,  the  user  inserts  a  blank  card  to  terminate  execution. 

2.3  SAMPLE  INPUT  SHEET  AND  OUTPUT 

Shown  below  is  a  typical  input  sheet  for  POST  showing  the  UNIVAC 
1108  "ADD"  cards  and  the  last  (blank)  card.  In  the  example  shown,  two  tests 
are  processed  U00100  and  U10105.  Following  the  sample  input  sheet  is  the 
actual  output  generated  by  POST. 
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2.4  The  listing  of  POST  appears 
Independent,  except  that  It  expects  the 
This  format  is  exactly  that  used  by  the 


below.  Note  that  It  is  machine 
Input  to  be  in  a  specific  format. 
NL001  code. 
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APPENDIX  B 


CHARACTERIZATION  CODE  NL002 


1.0  INTRODUCTION 

The  computer  code  NL002  is  a  nonlinear  thermoviscoelastlc  charac¬ 
terization  code  that  accepts  pressurized  uniaxial,  biaxial  and  shear  data 
from  code  NL001.  Dozens  of  experiments  of  this  type  having  arbitrary 
deformation-temperature  time  histories  can  be  simultaneously  fit  to  con¬ 
stitutive  equations  using  NL002  which  performs  both  a  distortlonal  and 
bulk  characterization  as  described  in  the  text  of  this  report. 

2.0  PROGRAM  INSTRUCTIONS 

The  purpose  of  this  section  is  to  present  the  basic  program  variables, 
their  definitions  and  use  ir.  the  NL002  code,  a  sample  problem  and  a  listing 
of  the  program  following  the  sample  problem  input  sheet. 

2.1  BASIC  VARIABLES 

A  Coefficient  in  e  A*DIL/EOCT  which  is  a  modulus  multiplier 

AT  Shift  function  -  input 

AVGAT  Interpolated  value  of  shift  function  at  a  specific  arbitrary  point 

p 

B  Coefficient  of  t  in  power  law  term  (1  +  Bt)  -  input 

BB  Solution  vector  in  regression  analysis 

BCAL  Calculated  value  of  dilatational  stress 

BETA  Thermal  volumetric  expansion  coefficient  -  input 

BO  Regression  coefficients  for  bulk  term  -  calculated 

BULK  Bulk  modulus  -  input 

CO  Regression  coefficients  for  shear  term  -  calculated 

D  Slope  of  shift  function  vs  temperature  curve  -  calculated 

DEV  Local  percent  error 


0-1 


DEVI 

Final  local  error  on 

DEV  2 

Final  local  error  on  o ^ 

DI 

Working  variable 

DIL 

Corrected  dilatation  AV/V  - 

0 

input 

DR1 

Calculated  strain  rate  E^  - 

input 

DR2 

Calculated  strain  rate  E ^  ~ 

input 

DR12 

Calculated  strain  rate  E^  - 

input 

DT 

Time  step  -  input 

DTR 

Reduced  time  step  -  calculated 

E0CT 

Octahedral  strain  -  input 

Ell,  E22, 
E33 

Principal  strains  -  input 

E12 

Shear  strain  -  input 

F  A  coefficient  defined  as  F  ■  e  [A*DIL/EQCT] 

E(VT 

G  A  coefficient  defined  as  G  »  (1  - 

IFF  F  coefficient  flag;  set  IFF(I)  -  1  to  multiply  Ith  term  in 

regression  by  the  coefficient  F 

IFG  G  coefficient  flag;  set  IFG(I)  ■  1  to  multiply  Ith  power  law 

term  by  the  coefficient  G 

IFN  Norm  selector;  selects  up  to  six  norm  terms,  10th,  20th,  30th, 

40th,  50th,  or  60th  depending  upon  whether  or  not  IFN(I)  is  1,  2, 
3,  4,  5,  or  6  respectively. 

INTI,  INT2,  Hereditary  integral  terms  -  calculated 

INT12 

INV1,  INV2,  First,  second,  and  third  strain  invariants  respectively  -  input 

INV3 


D-2 


IPG 


Page  councer 


K0DE  Test  type  indicator  -  input 

KTEMP  Constant  or  variable  temperature  test  indicator  -  input 

LAST  Input  data  delimiter 

LINE  Output  line  counter 

MATID  Material  identification  -  input 

NAME  Run  title  printed  at  top  of  each  page 

NAT  Number  of  values  of  shift  function  Vs  temperature 

NBULKS  Number  of  terms  in  series  for  bulk  term  regression 

presently  fixed  at  14 

NDP 

NORM  The  Lebesque  norms  of  Eq£T  -  input 

NORMF  Infinite  norm  of  "  input 

NORMS  Number  of  norm  terms  in  series  for  shear  term  regression  - 

limited  to  a  maximum  of  6 

P 

NPL  Number  of  power  law  terms  of  the  form  (HBt) 

NTERMS  Number  of  terms  in  series  for  shear  term  regression  -  limited 

to  a  maximum  of  14. 

NTEST  Test  number  in  this  run 

NT0T  Total  number  of  data  points  processed  this  run,  i.e.,  END? 

N1  First  data  point  used  in  regression  analysis 

p 

P  Exponent  in  power  law  term  (HBt) 

PRES  Test  pressure,  psig 

R  Primary  matrix  of  regression  analysis 


RATE 


Observed  strain  rate  -  input 


RES  ID 


RESID1 

RESID2 

SCAL 

STD 

STRESS 

STRN 

STRUE 

SUM 

SUM1 

SUM2 

SI,  S2 

S1CAL 

S2CAL 

T 

TEMP 

TITLE 

TR 

TSHFT 
XBAR 
XF,  XF1, 
XN 

Xt>,  XQ 
Z 


Interim  variable  used  In  calculating  atandard  deviations 

Calculated  value  of  stress  In  shear  regression 
Standard  deviation 
Observed  stress 

Interim  strain  storage  variable 
True  stress 

Interim  variables  used  in  calculating  average  error  and 
standard  deviation 

True  stress,  STRUE  -  S1-S2 

Calculated  values  of  SI  and  S2 

«« 

Time  In  minutes 
Test  temperature 
Test  Identification  -  input 
Reduced  time  -  calculated 

Temperature  at  which  shift  function  is  Input 
Average  error 

XF2  Terms  in  series  used  for  regression  analysis 

Power  to  which  norm  term  is  raised 
Working  variables 
Exponent  of  G  coefficient 


0-4 


2.2  PROGRAM  INPUT 

This  section  lists  the  required  Input  for  NL002  on  a  card  by 
card  basis.  The  format  for  each  card  Is  shown  In  parentheses. 

Card  1  (20A4)  NAME 

Card  2  (2E10.0)  A,  Z 

Card  3  (715)  IFG(I),  1-1,7 

Card  4  (1415)  IFF(I)  1  -  1,  14 

Card  5  (0615)  IFN(I) ,1  -  1,6 

Card  6  (6E10.0)  XN(Ie)  I  -  1,6 

Card  7  (15)  NPL 

Card(s)  (8E10)  (P(I)X  I  -  1  NPL) 

Card(s)  9  (2E10.0,  15)  TSMFT(I) ,  AT(I) ,  LAST.  Where 

LAST  *  0  or  blank  If  another  TSMFT  vs  AT  pair  follows,  and 
LAST  +  0  if  this  is  last  pair. 

Card(s)  10  "ADD"  cards  to  add  elements  from  NL001  output  data  tape,  as 
many  as  desired 

Card  11  A  blank  card  to  indicate  end  of  data 

The  above  cards  are  shown  on  the  coding  sheet  below. 

2.3  SAMPLE  PROBLEM 

Shown  below  is  an  input  sheet  for  a  simple  but  typical  problem 
using  13  tests  of  77°  data,  as  indicated  in  the  title,  Card  1.  Card  2  sets 
A  -  Z  »  1.0.  Card  3  indicates  that  *11  power  law  terms  are  multiplied  by  G. 
Card  4  indicates  that  all  terms  in  the  shear  term  series  are  to  be  multiplied 
by  F.  Card  5  says  to  take  6  norm  terms,  but  to  use  the  10th  norm  for  each 
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APPENDIX  E 


CHARACTERIZATION  CODE  NL003 

1.0  INTRODUCTION 

The  computer  code  NL003  Is  Identical  In  Input  and  operation  as 
NL002  and  performs  the  Identical  distortional  characterization.  Nl003 
differs  from  NL002  in  the  bulk  characterization.  In  this  code  the 
dilatation  is  expressed  In  terms  of  the  mean  pressure  and  the  octahedral 
shear  strain  as  discussed  In  the  text  of  this  report. 

Since  the  code  operation  is  identical  only  a  sample  problem  and 
the  listing  is  provided  in  this  appendix. 
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APPENDIX  F 


SUBPROGRAMS 

1.0  INTRODUCTION 

This  appendix  lists  the  additional  subprograms  called  for  In  code 
NL001,  NL002,  NL003,  POST,  and  LINVIS.  These  routines  are  used  for  paqe 
counting,  headings  and  matrix  inversion  and  are  common  to  these  codes. 
Since  these  subroutines  need  no  input  from  the  user  only  the  listings 
are  presented. 
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APPENDIX  6 

COMPUTERIZED  CHARACTERIZATION  PROCEDURES  FOR  LINEAR 
VISCOELASTIC  MATERIALS  USING  ARBITRARY  DEFORMATION  HISTORIES 


By 

Richard  J.  Farris 


1.0  INTRODUCTION 


The  purpose  of  this  appendix  Is  to  provide  an  Improved  means  of  linear 
viscoelastic  characterization.  The  computer  code  presented  In  this  appendix 
was  the  first  code  developed  on  this  contract  and  formed  the  basis  for  the 
nonlinear  characterization  codes.  The  code  calculates  the  best  fit  dlstor- 
tional  stress-strain  relation  for  mixed  uniaxial,  biaxial  and  shear  tests 
having  complex  deformation  histories  Including  transient  temperature  histories. 
The  representation  used  Is 


o(0)(t)  = 


t 

J'G(t'-c')e(D)(Odt 


where  =  dlstortlonal  stress 

=  distortional  strain 


Thermal  effects  are  Included  In  that  the  distortional  strain  depends 
upon  thermal  dilatation  and  the  reduced  time,  t ' -4 ' ,  depends  upon  temperature 
in  the  usual  manner 


t'-s’  = 


dx/Aj( t) 


The  kernel  function  G(t)  can  be  represented  as  either  a  Prony  series 
or  as  a  power  law  series  depending  upon  users  preference. 

In  the  Prony  series  representation  the  kernel  function  becomes 

M 

G(t)  =  A,  +  (3) 

1  1=2  1 


and  in  the  power  law  representation  the  kernel  function  is 


n  ^ 

G(t)  =  A]  +  1^(1  +  Bft)  ’ 


G-l  0J 


The  code  requires  as  Input  parameters  the  time -temperature  shift 
function  and  a  family  of  exponents  e-j ,  1  <  15,  for  the  Prony  series 
representation.  The  code  then  computes  tEe  best  fit  linear  coefficients 
Ai  based  on  a  relative  regression  analysis  as  discussed  In  the  text  of 
this  report.  Once  having  determined  the  coefficients  the  program  can  be 
used  to  calculate  the  response  to  any  history  of  Interest. 

There  are  many  advantages  to  this  type  of  characterization  method 
some  of  which  are  given  below. 

•  All  of  the  data  obtained  on  a  material  can  be  used  In  the 
characterization. 

•  The  method  can  provide  the  best  approximation  to  any  complex 
history. 

•  Provides  the  best  fit  to  all  of  the  data  not  just  a  relaxation 
test. 

•  Provides  direct  comparisons  between  observed  and  predicted  data. 

•  The  method  is  much  more  accurate  than  usual  methods  and  provides 
statistical  information  regarding  accuracy  and  variability. 

•  Eliminates  the  need  for  special  tests  and  test  equipment,  such 
as  stress-relaxometers,  since  characterization  can  be  carried 
out  from  any  experiment. 

•  It  provides  a  practical  method  for  linear  and  nonlinear  visco¬ 
elastic  characterization. 

As  an  example  of  the  application  of  the  method  the  following  figures 
are  included.  Figure  1  demonstrates  the  excellent  results  obtained  when  a 
five  term  Prony  series  Is  fit  to  a  ramp  loaded  stress-relaxation  experiment 
for  a  composite  propellant.  From  this  experiment,  one  might  get  the  mis¬ 
leading  impression  that  linear  viscoelasticity  was  a  good  approximation  to 
propellant  response.  Figure  2  shows  the  linear  viscoelastic  analytical 
continuation  (based  on  the  fit  from  Figure  1)  compared  to  the  experimental 
results  when  straining  is  again  commenced.  Obviously  the  material  Is  not 
linearly  viscoelastic.  The  remaining  figures  illustrate  the  excellent  fit 
obtained  in  Solethane  113,  an  unfilled  rubber,  on  a  series  of  experiments 
using  complex  histories  at  temperatures  from  -65°F  to  +150°F.  The  last 
figure  shows  a  comparison  of  predicted  transient  thermoviscoelastic  stress 
response  for  Solethane  113  with  experimental  results.  The  Solethane  charac¬ 
terization  was  carried  out  using  15  Prony  series  terms  and  the  predictions 
were  within  +  10*  for  one  standard  deviation  demonstrating  that  linear  visco¬ 
elasticity  is  an  excellent  approximation  to  the  behavior  of  this  unfilled 
polymer. 
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COMPARISON  OF  LINEAR  VISCOELASTIC  PREDICTION  AND  FIRST 
STRESS-RELAXATION  FOR  A  COMPOSITE  SOLID  PROPELLANT 


COMPARISON  OF  EXPERIMENTAL  AND  ANALYTICAL  EXTENSIONS  USING  CHARACTERIZATION 
FROM  FIRST  RELAXATION  FOR  THE  COMPOSITE  SOLID  PROPELLANT 
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Figure  2 


Time,  minutes 


COMPARISON  OF  LINEAR  VISCOELASTIC  PREDICTI 
AND  EXPERIMENTAL  DATA  FOR  SOLITKANE  113 
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Time,  minutes 


COMPARISON  OF  LINEAR  VISCOELASTIC  PREDICTIONS 
AND  EXPERIMENTAL  DATA  FOR  SOLITHANE  113 


PAR  ISON  OF  LINEAR  VISCOELASTIC  PREDICTI 
AND  EXPERIMENTAL  DATA  FOR  SOL  IT KANE  113 


COMPARISON  OF  LINEAR  VISCOELASTIC  PREDICTIONS 
AM'  EXPERIMENTAL  DATA  FOR  SOLITHANE  113 


Time,  minutes 


COMPARISON  OF  LINEAR  VISCOELASTIC  PREDICTIONS 
AND  EXPERIMENTAL  DATA  rOR  SOLITKANE  113 


Time,  minutes 


COMPARISON  OF  IINEAR  VISCOELASTIC  PREDICTIONS 
AND  EXPERIMENTAL  CATA  FOR  SOL  I THANE  113 


minutes 


COMPARISON  OF  LINEAR  VISCOELASTIC  PREDICTIONS 
AND  EXPERIMENTAL  DATA  FOR  SOL ITHAT  E  113 
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Figure  11 


mi rutes 


COMPARISON  OF  UNEAR  VISCOELASTIC  PREDICTIONS  ARC  EXPERIMENTAL  DATA 
FOR  THERMOVISCOELASTIC  EXPERIMENT  ON  SOLITKANE  113 


Tin,e,  hours 


I'KOCKAM  INSTRUCTIONS 


2.u 


This  section  presents  a  description  of  the  basic  varlabl  \s  of  the 
I.INVIS  linear  viscoelastic  characterization  code.  It  Includes  a  data  sheet 
showing  the  input  variables  and  formats,  a  sample  problem  and  a  listing  of 
the  program. 

2.  1  HAS IC  VARIABLES 


AT 

AVCAT 

B 

BB 

BETA 

CO 

1) 

nr.v 

DKVRAT 

l)T 

El,  E2 

E 

C 

12 

n 

[A 


Vector  of  shift  function  -  input 
Average  value  of  AT  -  calculated 


NN 


Constant  in  the  power  series  term  (1  +  Bt)  or  the  exponent 
in  the  prony  series  term  e  -  input 

Column  vector  in  regression  equation 

Volumetric  expansion  coefficient  -  input 

Regression  coefficients  -  calculated 

Slope  of  shift  function  -  calculated 

Percent  error  between  calculated  and  observed  stress 

Deviation  strain  rate  -  calculated 

Time  increment  -  input 

Principal  normal  strains  -  calculated 

Working  variable 

Inverse  matrix 

Set  to  1  to  include  a  linear  elastic  term  in  the  series  charac¬ 
terization  otherwise  leave  blank  or  zero  -  input 

Number  of  power  or  prony  series  terms  -  input 

12  +  13  -  calculated 


16 

INV1 

JL 

KFL 

KODE 

LLL 

NN 

NEXP 

NTEMP 

NTESTS 

RATE 

SIG1 

STD 

STRN1 

T 

TEM 

TEMP 

TIME 

TITLE 

TR 

TTEMP 

V 

VOLUME 


Equal  to  NTEMP  -  Calculated 

First  strain  Invariant  -  input 

Total  number  ot  input  data  points  -  calculated 

Number  of  input  data  points  in  a  particular  test  -  calculated 

Test  type:  •  0  fir  uniaxial;  **  1  for  shear;  -  2  for  biaxial 
stop  -  input 

Characterization  type:  »  0  for  power  series;  -  1  for  prony 
series  -  input 

NN 

Exponent  in  the  power  series  term  (1  +  Bt)  -  input 
Number  of  data  points  for  this  test  -  input 

ilumber  of  shift  function  vs  temperature  pairs  -  input,  20  max. 

Number  of  tests  for  this  run  -  input 

Incremental  strain  rate  -  input 

Working  array  to  store  all  stress  input 

Standard  deviation  -  calculated 

Same  as  SIG1  except  for  strain  -  calculated 

Time  -  calculated 

Same  as  SIG1  except  for  temperature  -  calculated 
Temperature  at  a  data  point  -  input 
Output  time  -  calculated 

80  column  alphanumeric  identification  -  input 
Reduced  time  -  calculated 
Shift  function  input  temperature  -  input 
Working  matrix 

Same  as  SIG  )  except  for  INV1  -  calculated 
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X 


Dummy  array 


XBAR  Average  error  -  calculated 

XF  Union  of  the  X  sets 

Y  Observed  stress  -  input 

YCAL  Calculated  stress  -  output 

2 . 2  DATA  INPUT 

The  input  variable  required  by  LINVIS  are  shown  below,  card  by 
card,  with  the  format  shown  in  parentheses. 

Card  1  (20A4)  TITLE 
Card  2  (E10.0)  BETA 

Card  3  (515)  12,  13,  LLL,  NTESTS ,  NTEMP 
Card(s)  4  (2E10.0)  NN,  8 
Card(s)  5  (2E10.0)  TTEMP ,  AT 
Card  6  (215)  NEXP,  K0DF. 

Card(s)  7  (5E10.0)  TEMP,  RATE,  DT,  Y,  INV1 

Cards  1  through  5  are  read  only  once  per  run. 

Cards  6  and  7  are  repeated  for  each  test. 

The  above  cards  are  shown  on  the  data  sheet  below. 

2.3  Shown  below  are  the  sample  input  sheets  for  a  typical  problem. 

By  referring  to  the  master  data  sheet,  above,  the  input  is  self  explanatory. 
The  output  sheets  for  this  sample  problem  follow  the  data  sheets. 

2.4  PROGRAM  LISTING 

A  complete  source  listing  follows  the  sample  problem  output. 
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APPENDIX  H 

SAMPLES  OF  EXPERIMENTAL  DATA 


This  appendix  lists  some  of  the  experimental  data  obtained  on 
ANB-3335-1  propellant  at  0°F.  These  data  are  from  "pressurized"  uniaxial 
and  biaxial  experiments.  In  total  over  250  experiments  of  the  type  Illus¬ 
trated  were  performed  on  this  propellant  at  temperatures  from  -65°F  to 
+150°F  and  pressures  from  atmospheric  to  1000  psl.  After  processing  the 
raw  data  through  the  Preprocessor  Code  NL001,  calculated  stress-strain- 
invariant  histories  are  stored  on  magnetic  tape.  The  experimental  data 
listed  In  this  appendix  represents  only  a  partial  listing  of  the  calculated 
Information  from  each  experiment.  These  listings  were  obtained  using  the 
Post  Processor  Code. 

These  data  are  Illustrated  In  this  appendix  to  demonstrate  the  type 
of  experimental  data  obtained  and  used  on  this  contract.  Arrangements  can 
be  made  to  obtain  a  copy  of  the  complete  magnetic  tape  data  file  for  those 
who  are  Interested  In  testing  theory  against  real  data. 
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